AD-AU9  704 
UNCLASSIFIED 


MASSACHUSETTS  INST  OF  TECH  CAMBRIDGE  CENTER  FOR  MATE— ETC  F/6  11/6  V 

MOLECULAR-ORBITAL  BASIS  FOR  SUPERCONDUCTIVITY  IN  HIGH-  AND  LOW— ETC <U) 

SEP  82  K  H  JOHNSON#  R  P  MESSMER  N00014-B1-K-0499 

TR-4  NL 


Tf< 

00 

i> 

Cft 


Unclassified 

/  ■  —  .-■■■  - 
*CCuftlTY  CLA»1FICATiQh  OF  THIS  PAGE  fWft»o  Pete  Entered) 

I  REPORT  DOCUMENTATION  PAGE 


REPORT  mUVICK 


1 4.  title  (tntt  Submit) 


|2.  GOVT  ACCESSION  NOJ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 
RECIPIENT’S  CATALOG  NUMGER 


S.  TYPE  OP  REPORT  t  PERIOO  COVERED 


Molecular-Orbital  Basis  for  Superconductivity 
In  High-  and  Low-Dimensional  Metals 


Interim 


17.  AuTmOR<*J 


«.  PERFORMING  OAG.  REPORT  NUMGER 


I.  CONTRACT  OR  GRANT  NUMGERft) 


K.  H.  Johnson  and  R.  P.  Messmer  N00014-81-K-0499 

>.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS  10.  JJgORAM  ELEMENT  PROJECT.  T 

f  Ant  4  v  U  N I  T  NUMItA# 

Center  for  Materials  Science  and  Engineering, 

M.I.T.,  Cambridge,  Massachusetts  02139  Task  No.  Nr  056-757 

II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS  '*•  REPORT  OATE 

Office  of  Naval  Research  September  24,  1982 

Department  of  the  Navy  11  NUMBC*  or  PA<?«fs 

_ Arlington.  Virginia  22217 _ _ _ 

14.  MONITORING  AGENCY  NAME  «  ADDRESS (U  dllltttnl  from  Controlling  Olllct)  IS.  SECURITY  CLASS,  (ol  ihlt  import) 


ID.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  t  WORK  UNIT  NUMGERS 


Task  No.  Nr  056-757 

12.  REPORT  OATE 

September  24,  1982 

13.  NUMBER  OF  PAGES 

98 


15«.  oeclassification/Wwngrading 
SCHEDULE 


Ms.  distribution  statement  (of  thia  Report) 


kELECTE) 
„  OCT  1  1982 


Approval  for  public  release; 
distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  thm  abetract  entered  in  Block  20,  if  difforont  from  Roport) 


j  is.  SUPPLEMENTARY  notes 


13.  KEY  WORDS  (Continue  on  reveree  elde  If  neceeeary  and  Identify  by  block  number) 

Molecular  Orbital 
Superconductivity 

)  High  and  Low  Dimensional  Metals 


20.  VflfeTRACT  (Continue  on  ravoraa  aide  If  nacaaaary  and  Identify  by  block  number) 

^ real-space  molecular-orbital  description  of  electronic  wave  functions 
which  are  postulated  to  be  the  precursors  of  the  superconducting  state  In 
high-  and  low-dimensional  metals  is  presented,  baspd  on  self-consistent 
X-alpha  scattered-wave  (SCF-Xa-SW)  molecular-orbital  calculations  for 
clusters  representing  the  local  molecular  environments  In  these  materials. 
It  is  shown  that  there  is  a  persistent  correlation  between  the  occurrence 
of  superconductivity  in  a  material  and  th4  existence  of  spatially  y 


DD  1473  eoition  of  i  nov •»  is oosolete  Unclassified 

J'lQll  02-LQ31J-6601  Q  Q  Q  SECURITY  CLASSIFICATION  OF  TH|»  PAGE  Dbtb  Tntrnrt*) 


_ Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PACE  <*hm  Dal*  Enl*r«Q 


^delocalized  molecular  orbitals  at  the  Fermi  energy  which  are  bonding  within 
and  anti  bonding  between  “layers"  or  "tubes"  of  overlapping  atomic  orbitals  that 
span  many  atoms,  forming  a  type  of  "electron  network"  at  the  Fermi  energy,  as 
exemplified  by  pir  "layered"  molecular-orbital  topologies  In  A1  and  (TMTSF)2PFfi, 
and  by  d6  "tubular"  molecular-orbital  topologies  in  Nb  and  Nb3Sn.  This  c  D 
description  of  the  precursor  superconducting  state  is  consistent  with  the 
original  conjectures  of  London  that  the  superconducting-state  wave  function  Is 
"molecular"  in  nature,  "rigid"  In  character,  and  of  wide  spatial  extent,  from 
which  observed  physical  properties  (e.g.,  diamagnetism  and  nondlsslpatlve 
electrical  currents)  of  the  superconducting  state  logically  follow.  The 
molecular-orbital  model  is  further  shown  to  be  consistent  with  Coopfrjs  con¬ 
cept  of  electron  pairing  in  the  superconducting  state  through  a  net 'attractive 
el ectron-el ectron  interaction  but  differs  from  the  Bardeen-Cooper-Schrreffer 
(BCS)  theory  of  superconductivity  in  attributing  the  pairing  to  "valence-bond- 
like"  electron  occupation  of  the  layered  or  tubular  molecular  orbitals  at  the 
Fermi  energy,  coupled  with  lattice  ion  displacements  through  a  dynamic  Jahn- 
Teller  effect,  rather  than  to  the  absorption  and  emission  of  virtual  phonons. 
However,  recent  simple  physical  arguments  presented  by  Weisskopf  for  the 
formation  of  Cooper  pairs  suggest  that  the  two  points  of  view  are  not  in¬ 
compatible.  Approximate  formulae  for  calculating  the  electron  pair  binding 
energy  and  superconducting  transition  temperature,  analogous  to  the  BCS 
formulae  but  expressed  entirely  in  terms  of  two  simple  molecular-orbital 
parameters,  are  derived  and  applied  to  various  superconductors.  Moreover, 
through  an  elementary  phase-space  argument,  the  molecular-orbital  description 
of  Cooper  pairs  is  shown  to  be  qualitatively  consistent  with  Schafroth's 
description  of  electron  pairs  as  "quasi -molecules"  undergoing  Bose-Einstein- 
like  condensation  to  the  superconducting  state.  The  above  molecular-orbital 
criteria  for  superconductivity  are  diametrical  to  those  for  the  occurrence  of 
local  magnetic  moments  and  ferromagnetism,  namely,  the  existence  of  spatially 
localized,  exclusively  anti  bonding  (e.g.,  do*,  dir*,  or  d6*)  molecular  orbitals 
at  the  Fermi  energy.  Furthermore,  in  contrast  to  the  ordering  of  electron 
spins  on  atomic  sublattices  in  conventional  anti  ferromagnetism,  the  valence- 
bond-like  correlation  of  Cooper-Schafroth  electron  pairs  at  the  Fermi  energy 
in  composite  bonding-anti  bonding  layered  or  tubular  molecular-orbital 
components,  spatially  delocalized  between  the  atoms,  corresponds  to  a  type  of 
conduction-electron  "anti ferromagnetism"  or  spin-density  wave,  thus  offering 
an  explanation  for  the  occurrence  of  such  anti ferromagnetism  in  super¬ 
conductors  such  as  (TMTSF^PFc  and  NbjSn.  These  molecular-orbital  criteria 
therefore  provide  a  conceptual  basis  for  understanding  the  generally  mutually 
exclusive  incidence  of  superconductivity  and  magnetism  among  the  elements  of 
the  periodic  table,  although  they  can  also  be  used  to  explain  the  occasional 
coexistence  of  superconductivity  and  ferromagnetism  or  anti ferromagnetism  in 
some  materials.  The  molecular-orbital  model  for  superconductivity  complements 
BCS  theory  in  that  it  permits  the  prediction  of  which  materials  are  likely  to 
be  superconductors  and  which  are  not  entirely  on  the  basis  of  the  molecular- 
orbital  topology  at  the  Fermi  energy.  It  offers  an  explanation  of  the  vanish¬ 
ing  isotope  effect  in  certain  superconductors.  It  contributes  to  the 
clarification  of  the  issue  of  non-phonon  mechanisms  of  superconductivity.  It 
provides  insight  into  the  correlations  of  superconductivity  with  other  physical 
properties,  such  as  lattice  instabilities.  It  is  readily  applicable  to  super¬ 
conductors  lacking  long-range  crystalline  order,  suph  as  amorphous  alloys  and 
small  particles.  Finally,  the  molecular-orbital  approach  can  be  used  to 
explain  in  simple  terms  why  some  materials  (e.g.,  Cu,  Ag,  and  Au)  are  neither 
superconducting  nor  magnetic,  why  certain  quasi -one-dimensional  organic  solids. 
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SUMMARY 

A  real-space  molecular-orbital  description  of  electronic  wave  functions  which 
are  postulated  to  be  the  precursors  of  the  superconducting  state  in  high-  and  low¬ 
dimensional  metals  is  presented,  based  on  self-consistent  X-alpha  scattered-wave 
(SCF-Xa-SW)  molecular-orbita1  calculations  for  clusters  representing  the  local 
molecular  environments  in  these  materials.  It  is  shown  that  there  is  a  persistent 
correlation  between  the  occurrence  of  superconductivity  in  a  material  and  the 
existence  of  spatially  delocalized  molecular  orbitals  at  the  Fermi  energy  which 
are  bonding  with  !n  and  anti  bonding  between  "layers"  or  "tubes"  of  overlapping 
atomic  orbitals  that  span  many  atoms,  forming  a  type  of  "electron  network"  at 
the  Fermi  energy,  as  exemplified  by  pm  "layered"  molecular-orbital  topologies 
in  A1  and  (TMTSF^PF^,  and  by  do  "tubular"  molecular-orbital  topologies  in 
Nb  and  Nb3Sn.  This  description  of  the  precursor  superconducting  state  is 
consistent  with  the  original  conjectures  of  London  that  the  superconducting- 
state  wave  function  is  "molecular"  in  nature,  "rigid"  in  character,  and  of 
wide  spatial  extent,  from  which  observed  physical  properties  (e.g.,  diamagnetism  and 
nondissipative  electrical  currents)  of  the  superconducting  state  logically  follow. 
The  molecular-orbital  model  is  further  shown  to  be  consistent  with  Cooper's 
concept  of  electron  pairing  in  the  superconducting  state  through  a  net  attractive 
elegtron-electron  interaction  but  differs  from  the  Bardeen-Cooper-Schrieffer 
(BCS)  theory  of  superconductivity  in  attributing  the  pairing  to  "valence-bond- 
like"  electron  occupation  of  the  layered  or  tubular 
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molecular  orbitals  at  the  Fermi  energy,  coupled  with  lattice  Ion  displacements 
through  a  dynamic  Jahn-Teller  effect,  rather  than  to  the  absorption  and  emission 
of  virtual  phonons.  However,  recent  simple  physical  arguments  presented  by 
Weisskopf  for  the  formation  of  Cooper  pairs  suggest  that  the  two  points  of  view 
are  not  incompatible.  Approximate  formulae  for  calculating  the  electron 
pair  binding  energy  and  superconducting  transition  temperature,  analogous  to  the 
BCS  formulae  but  expressed  entirely  in  terms  of  two  simple  molecular-orbital 
parameters,  are  derived  and  applied  to  various  superconductors .  Moreover, 
through  an  elementary  phase-space  argument,  the  molecular-orDital  description 
of  Cooper  pairs  is  shown  to  be  qualitatively  consistent  with  Schafroth's 
description  of  electron  pairs  as  "quasi-molecules"  undergoing  Bose-Einstein- 
like  condensation  to  the  superconducting  state.  The  above  molecular- 
orbital  criteria  for  superconductivity  are  diametrical  to  those  for  the 
occurrence  of  local  magnetic  moments  and  ferromagnetism,  namely,  the  existence 
of  spatially  localized,  exclusively  anti  bonding  (e.g.,  dc*,  d-*,  or  d5*) 
molecular  orbitals  at  the  Fermi  energy.  Furthermore,  in  contrast  to  the  order¬ 
ing  of  electron  spins  on  atomic  sublattices  in  conventional  anti  ferromagnetism, 
the  valence-bond-like  correlation  of  Cooper-Schafroth  electron  pairs  at  the 
Fermi  energy  in  compos i te  bonding-antibonding  layered  or  tubular  molecular- 
orbital  components,  spatially  delocalized  between  the  atoms,  corresponds  to 
a  type  of  conduction-electron  "anti ferromagnetism"  or  spin-density  wave,  thus 
offering  an  explanation  for  the  occurrence  of  such  anti  ferromagnetism  in  super¬ 
conductors  such  as  (TMTSF^PFg  and  Nb^Sn.  These  molecular-orbital  criteria 
therefore  provide  a  conceptual  basis  for  understanding  the  generally  mutually 
exclusive  incidence  of  superconductivity  and  magnetism  among  the  elements  of 
the  periodic  table,  although  they  can  also  be  used  to  explain  the  occasional 
coexistence  of  superconductivity  and  ferromagnetism  or  anti  ferromagnetism  in  some 
materials.  The  molecular-orbital  model  for  superconductivity  complements 
BCS  theory  in  that  it  permits  the  prediction  of  which  materials  are  likely 


to  be  superconductors  and  which  are  not  entirely  on  the  basis  of  the  molecular- 
orbital  topology  at  the  Fermi  energy,  it  offers  an  explanation  of  the  vanish¬ 
ing  Isotope  effect  in  certain  superconductors.  It  contributes  to  the 
clarification  of  the  issue  of  non-phonon  mechanisms  of  superconductivity.  It 
provides  insight  into  the  correlations  of  superconductivity  with  other  physical 
properties,  such  as  lattice  instabilities.  It  is  readily  applicable  to  super¬ 
conductors  lacking  long-range  crystalline  order,  such  as  amorphous  alloys 
and  small  particles.  Finally,  the  molecular-orbital  approach  can  be  used 
to  explain  in  simple  terms  why  some  materials  (e.g.,  Cu,  Ag,  and  Au)  are 
neither  superconducting  nor  magnetic,  why  certain  quasi -one-dimensional  organic 
solids,  such  as  TTF-TCNQ,  are  not  superconductors  while  others,  such  as 
(TMTSF)2PF6,  are  superconductors,  and  to  assess  critically  the  likelihood  of 
superconductivity  in  certain  other  types  of  substances  (e.g.,  other  types  of 
organic  solids;  metallic  hydrogen  at  attainable  high  pressures),  and  to 
suggest  ways  of  systematically  improving  existing  classes  or  synthesizing  novel 
classes  of  superconducting  materials. 

I .  INTRODUCTION 

The  Bardeen-Cooper-Schrieffer  (BCS)  theory  of  superconductivity  [1]  ascribes 

the  onset  of  the  superconducting  state  of  a  crystal  at  the  transition  temperature, 

Tc>  to  electrons  attractively  paired  via  virtual  phonons  of  the  lattice.  The 

BCS  formula  for  T  is 
c 

Tc  n.  0Q  exp  [-l/N(0)V]f  (1) 

where  6q  is  the  lattice  Debye  temperature,  V  is  the  attractive  potential  between 
electrons  in  Cooper  pairs  [2]  induced  by  the  electron-phonon  interactions,  and 
N( 0)  is  the  electronic  density  of  states  for  one  spin  at  the  Fermi  energy. 

The  BCS  theory  together  with  its  strong-coupling  extension  has  been  eminently 
successful  in  accounting  for  the  physical  properties,  e.g.,  nondissipative 
current,  diamagnetism,  and  thermodynamics  of  the  superconducting  state,  and 


has  correlated  many  expermental  data  In  terms  of  a  few  basic  parameters. 
Nevertheless,  it  has  often  been  emphasized  by  Matthias  [3]  that  the  BCS  theory 
and  formal  extensions  thereof  do  not  satisfactorily  explain  the  observed 
dependence  of  superconductivity  on  crystal  structure  and  chemistry  (especially 
for  transition  metals,  alloys,  and  compounds)  and  are  not  very  useful  for 
predicting  which  materials  should  be  superconducting  and  which  should  not. 

The  existence  of  superconducting  materials  possessing  only  short-range  structural 
order,  such  as  superconducting  amorphous  alloys  [4]  and  superconducting  metal 

O 

particles  down  to  50A  in  size,  [5,6]  would  appear  to  be  more  conveniently 
viewed  from  a  local  "real -space"  chemical  approach,  than  by  traditional  concepts 
of  long-range  cyrstalline  order  and  momentum  (ic)  space,  which  lead  to  coherence 
lengths  of  the  superconducting  state  that  usually  exceed  the  short-range  order 
and  electron  mean  free  path  characteristic  of  amorphous  and  small -particle 
superconductors.  This  emphasizes  the  desirability  of  having  a  local  chemical¬ 
bonding  or  real-space  molecular  description  of  the  superconducting  state  inorder 
to  complement  BCS  theory.  Indeed,  London  [7,8]  in  his  phenomenological  approach 
to  superconductivity  discusses  the  possibility  of  developing  a  molecular 
description  of  the  superconducting  state  'see  Chapter  E  of  Ref.  8),  and  Slater  [9] 
in  an  early  attempt  at  describing  superconductivity  discusses  the  nature  of 
the  spatial  character  of  the  superconducting-state  wave  function.  With 
speculations  that  mechanisms  other  than  electron-phonon  coupling  can  attractively 
pair  electrons  in  the  superconducting  state,  [3,10,11]  a  molecular  criterion 
that  accounts  for  the  known  chemical  trends  in  the  occurrence  of  superconductivity 
would  be  a  useful,  tool  in  the  continuing  effort  to  identify  new  classes  of 
superconducting  materials. 

The  systematic  trends  of  T£  with  valence-electron-per-atom  ratio  [3]  and 
the  difficulty  of  calculating  BCS  parameters  for  transition  metals,  their 
alloys,  and  compounds  have  prompted  the  development  of  rudimentary  chemical 
models,  [12-14]  based  on  atomic  d  orbitals  and  afmed  primarily  at  the  estimation 
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of  Tc-  Unfortunately,  these  models  are  not  easily  generalized  to  non- 
transition-metal  superconductors  and  do  not  yield  a  real-space  molecular 
description  of  the  superconducting  state.  The  most  general  chemical  approach 
to  superconductivity  is  that  of  Krebs,  [15]  who  has  suggested  that  superconduc¬ 
tivity  is  possible  only  if  the  normal  chemical  bonding  system  in  the  crystal 
or  parts  of  the  crystal  permits  the  construction  of  a  molecular  wave  function 
which  for  at  least  one  space  direction  is  not  intersected  by  plane  or  conical 
nodal  surfaces,  and  if  the  corresponding  electron  band  is  not  fully  occupied. 
However,  Krebs  attempts  only  to  account  for  the  incidence  of  superconductivity 
through  qualitative  chemical  bonding  arguments  and  not  the  evaluation  of  a 
specific  electron-pairing  mechanism  and  T  .  Local  reformulations  of  BCS-type 
superconductivity  theory  for  crystals  in  terms  of  Wannier  functions  [16-18] 
should  also  be  mentioned,  although  these  approaches  do  not  establish  a  direct 
relationship  between  superconducti vity  and  the  nature  of  the  chemical  bond.  In 
the  latter  regard,  one  should  not  fail  to  ci  .e  the  paper  of  Pauling  [19]  in 
which  a  relationship  between  "resonating  valence  bonds"  in  metals  and  the  BCS 
electron-phonon  mechanism  for  superconductivity  is  proposed. 

In  this  paper,  we  present  a  real -space  molecular-orbital  description  of 
the  wave  functions  which  we  believe  to  be  the  precursors  of  the  superconducting 
state,  as  well  as  providing  simple  criteria  for  determining  the  occurrence  of 
superconductivity.  This  theoretical  model  for  superconductivity  makes  direct 
use  of  the  results  of  first-principles  self-consistent-field  X-alpha  scattered- 
wave  (SCF-Xa-SW)  molecular-orbital  calculations  for  clusters  representing  the 
local  molecular  environments  in  these  materials.  The  SCF-Xa-SW  method  [20]  is 
based  on  the  combined  use  of  the  Xa  density-functional  approximation  [21]  to 
electron-electron  exchange  and  correlation  and  the  multiple-scattered-wave 
technique  [22]  of  solving  the  SchrSdinger  self-consistent-field  molecular-orbital 
equations.  In  conjunction  with  the  "transition-state"  procedure  [20,21]  for 
calculating  electronic  transitions  between  orbitals,  this  method  leads  to  a 


description  of  both  the  ground  and  excited  electronic  states  of  a  system.  The 
calculation  of  different  molecular  orbitals  for  different  spins  allows 
one  to  consider  the  effects  of  magnetic  spin  polarization.  Along  with 
applications  to  a  wide  variety  of  polyatomic  molecules,  [23,24]  including 
biological  macromolecules  [25-27]  and  polymers.  [28]  this  approach  has  been  used 
successfully  to  elucidate  the  local  bulk  and  surface  electronic  structures 
and  properties  of  solids  [29-31]  e.g.,  crystalline  semiconductors  and  localized 
defects  and  impurities  therein,  [32-35]  amorphous  semiconductors,  [36]  crystalline 
and  amorphous  metals  and  alloys  [37-42]  including  magnetic  properties  [39-41] 
molecular  crystals,  [28]  and  chemisorption  on  metal  surfaces  [43,44]  along  with 
the  calculation  of  photoelectron  emission  intensities,  [45,46]  electron-scattering 
cross  sections,  [47]  and  "shake-up"  spectra  [48]  for  chemisorbed  species.  These 
theoretical  studies  collectively  suggest  that  short-range  order  and  local  chemical 
bonding  largely  determine  the  measurable  electronic  structures  and  properties  of 
many  types  of  materials,  whether  ordered  or  disordered. 

In  view  of  the  apparent  complementary  natures  of  superconductivity  and 
magnetism,  it  should  be  noted  that  spin-unrestricted  SCF-Xa-SVJ  cluster  molecular- 
orbital  studies  [39]  indicate  that  the  local  magnetic  moments  of  Fe  and  Mn 
impurities  in  crystalline  copper  above  the  Kondo  temperature  are  associated 
with  the  existence  of  exclusively  antibonding  Fe(d)-Cu(s,d)  molecular  orbitals  at 
the  Fermi  energy  (see  Figs.  10,  12,  and  13  of  Ref.  39).  Similar  spin-polarized 
cluster  models  for  ferromagnetic  bcc  crystalline  ct-iron  [41,30]  suggest  that  the 
onset  of  collective  magnetism  below  the  Curie  temperature  and  the  persistence  of 
"spin  clusters"  above  the  Curie  temperature  are  associated  with  the  existence 
of  exclusively  antibonding  Fe(d)-Fe(d)  molecular  orbitals  at  the  Fermi  energy 
in  the  non-spin-polarized  limit  (see  Fig.  14  of  Ref.  41).  Recent  spin-polarized 
SCF-Xa-SW  cluster  models  for  amorphous  iron-boron  and  iron-phosphorus  alloys  [42] 
are  also  in  good  quantitative  agreement  with  measured  magnetic  properties  and 
clearly  show  the  key  role  of  anti  bonding,  molecular  orbitals  at  the  Fermi  energy 


in  determining  how  alloying  alters  the  magnetic  properties  of  Iron.  For  the 
sake  of  completeness.  It  may  also  be  noted  that  the  SCF-Xa-SW  method  has 
been  applied  to  iron-containing  biological  molecules  such  as  ferredoxin  [26] 
and  hemoglobin  [27]  and  the  results  used  to  explain  observed  magnetic  properties 
(including  the  calculation  of  magnetic  hyperfine  parameters  [27])  as ‘functions 
of  temperature  and  oxygenation,  respectively.  Finally,  the  spin-unrestricted 
SCF-Xa-SW  method,  in  conjunction  with  the  transition-state  procedure,  [20,21] 
has  recently  been  employed  with  good  results  to  calculate,  from  first  principles, 
ferromagnetic  and  anti  ferromagnetic  exchange  coupling  constants  (Heisenberg 
exchange  integrals)  for  transition-metal  cluster  coordination  o  * :xes .  [49] 

The  most  general  finding  of  the  theoretical  studies  of  magnetic  stances 
summarized  above  is  that  for  the  local  and  collective  forms  of  ^tism 

investigated,  one  finds  the  occurrence  of  antibonding  molecular  cal s  near 
the  Fermi  energy  in  the  non-spin-polarized  limit  (see  further  discussion  in 
Section  VI). 

In  view  of  the  above  work,  it  was  thought  that  SCF-Xa-SW  cluster  molecular- 
orbital  studies  of  superconducting  materials  might  possibly  also  provide  some 
new  insight  into  the  nature  of  chemical  bonding  associated  with  the  super¬ 
conducting  state.  As  will  be  discussed  in  Section  IV,  it  appears  that  the 
molecular-orbital  criteria  for  the  occurrence  of  superconductivity  and  certain 
types  of  magnetism  are  diametrical,  thereby  offering  a  simple  explanation 
for  the  mutually  exclusive  occurrence  of  these  phenomena  in  a  material  at 
the  same  temperature  and  pressure.  The  orbital  criteria  for  superconductivity 
will  further  be  shown  to  be  consistent  with  Cooper’s  [2]  concept  of  electron 
pairing  through  a  net  attractive  electron-electron  interaction  (V  in  Eq.  1) 
but  differs  somewhat  from  the  generally  accepted  mechanism  for  superconductivity 
in  attributing  this  attractive  interaction  directly  to  real -space  chemical - 
bond-like  formation  at  the  Fermi  energy,  coupled  with  the  lattice  through  a 
dynamic  Jahn-Teller  effect,  rather  than  to  electron-phonon  coupling  per  se. 
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However,  arguments  will  be  presented  to  show  that  these  two  points  of  view 
are  not  necessarily  incompatible.  The  molecular-orbital  approach  presented 
herein  will  also  be  used  to  explain  simply  why  some  metals  (e.g.  Cu,  Ag,  and  Au) 
are  neither  superconducting  nor  magnetic.  Finally,  it  is  suggested  that 
molecular-orbital  criteria  might  provide  a  useful  theoretical  tool  for  ' 
predicting  new  superconducting  materials. 

II.  CLUSTER  MOLECULAR  ORBITALS  AND  BAND  STRUCTURE 

The  relationships  of  SCF-Xa-SW  cluster  molecular-orbital  energy  levels 
and  densities  of  states  to  the  corresponding  crystal  band  structures  and 
densities  of  states  calculated  from  Bloch's  theorem  and  measured  spectroscopically 
have  been  discussed  in  previous  publications  and  reports  for  the  examples  of 
copper,  [37,39]  nickel,  [37]  aluminum,  [38,43]  and  iron.  [41]  Nevertheless, 
in  order  to  clarify  the  use  of  cluster  molecular  orbitals  as  a  basis  for 
superconductivity,  we  point  out  some  specific  similarities  and  differences 
between  cluster  electronic  structure  and  crystal  band  structure. 

For  simplicity,  consider  first  the  example  of  an  isolated  benzene 
molecule,  the  p-n  orbitals  of  which  are  shown  schematically  on  the  le  t  *ide 
of  Fig.  1.  The  orbital,  which  is  occupied  by  two  spin-paired  electrons 

and  has  the  energy  level  ,  is  the  most  strongly  bonding  of  the  pn  molecular 
orbitals,  whereas  the  degenerate  ^  and  ^  orbitals,  which  together  are 
occupied  by  four  spin-paired  electrons  and  have  the  energy  level  z^  are 
relatively  weakly  bonding.  The  unoccupied  orbital  is  strongly  antibonding, 
while  the  unoccupied  rr^*  and  orbitals  are  weakly  antibonding.  Because 
of  the  one-dimensional  periodicity  of  the  carbon  atoms  around  the  benzene 
ring,  the  molecular  orbitals  can  be  expressed  in  symmetrized  Bloch  form,  [50] 
resulting  in  the  discrete  one-dimensional  "k-space"  band  profile  shown 
schematically  on  the  right  side  of  Fig.  1.  The  "Fermi  energy"  z^  separates 
the  occupied  states  from  the  unoccupied  ones.  In  this  particular  example. 
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there  is  a  simple  one-to-one  correspondence  between  the  "real  space" 
molecular  orbitals  and  "k-space"  band  structure,  as  indicated  by  the  connecting 
lines  in  Fig.  1.  It  has  often  been  noted  that  the  diamagnetic  "ring  currents" 
induced  in  benzene  and  other  aromatic  molecules  through  the  response  of  the 
piT  electrons  to  an  external  magnetic  field  are  nondisslpative  currents  similar 
in  some  (but  not  all)  respects  to  the  persfstent  currents  of  a  superconducting 
ring.  [8,51,52] 

On  the  other  hand,  the  molecular  orbitals  of  a  finite  cluster  representing 
the  local  molecular  environment  in  a  crystal  transform  according  to  the  irreducible 
representations  of  the  cluster  point  symmetry  group  and  therefore  are  not  generally 
equivalent  on  a  one-to-one  basis  to  Bloch  band  eigenstates  of  the  crystal, 
except  at  the  points  in  £-space  which  have  the  full  cluster  point  group  symmetry. 
For  example,  the  T(k=0)  Bloch  state  has  the  full  0^  point-group  symmetry  in 
a  cubic  crystal.  Nevertheless,  it  is  possible  to  construct  combinations  of 
degenerate  Bloch  states  of  different  k  (the  so-called  "star"  of  the  wave  vector) 
that  transform  according  to  the  irreducible  representations  of  the  cluster  point 
group.  [53]  From  this,  one  can  conclude  that  there  should  be  a  reasonably  close 
correspondence  between  cluster  and  crystal  densities  of  states.  Examples  of 
such  correspondence  will  be  found  in  Refs.  37,38,41,  and  43. 

As  an  illustrative  application  of  the  above  concepts  to  a  superconducting 
metal,  consider  the  example  of  aluminum.  In  Fig.  2,  the  SCF-Xa-SW  molecular- 
orbital  energy  levels  of  a  43-atom  cluster,  AlAl-^AIgA^,  representing  the 
local  environment  of  fee  crystalline  aluminum  up  to  third-nearest  neighbors  [38] 
are  compared  with  the  energy  bands  of  the  crystal  [54]  shown  along  high  symmetry 
directions  of  the  Brillouin  zone.  The  principal  nearest-neighbor  chemical 
bonding  components  (o,  -t,  etc.)  of  the  A1  s  and  p  orbitals  in  the  cluster 
are  indicated.  The  distribution  of  cluster  energy  levels  is  also  compared 
with  the  X-ray  photoelectron  spectrum  (XPS)  of  pure  crystalline  aluminum  in 
Fig.  3.  Good  agreement  between  the  cluster  results  and  the  XPS  spectrum  is 


obtained  with  respect  to  both  density  and  total  s,p  band  width  of  the 
cluster  energy  levels,  consistent  with  the  comparison  of  aluminum  cluster 
densities  of  states  and  ultraviolet  photoemission  spectra  (UPS)  described  in 
Ref.  38.  There  is  a  close  correspondence  between  individual  cluster  molecular 
orbitals  and  certain  81och  energy-band  eigenstates  or  combinations  thereof. 

For  example,  the  lowest  occupied  so-bonded ‘cluster  energy  level,  la.^,  shown 
in  Fig.  2  is  the  discrete  molecular-orbital  analogue  of  the  r1  (1c* 0)  s-band 
Bloch  state.  The  manifold  of  s,p  o-bonded  cluster  levels  in  the  energy  range 
-1.3  Ry  to  -0.7  Ry  corresponds  to  combinations  of  A-j  and  E-j  Bloch  states  in 
the  first  band  of  aluminum,  with  the  high  density  of  s,p  o-bonded  cluster 
levels  between  -0.85  Ry  and  -0.65  Ry  responsible  for  a  broad  maximum  in  the 
XPS  density  of  states  evident  in  Fig.  3.  The  high  density  of  po-  and  pir- 
bonded  cluster  levels  between  -0.65  Ry  and  the  Fermi  energy  is  associated 
with  combinations  of  ,  E^,  and  Z^  Bloch  states  in  the  second  band  of 
aluminum  and  is  responsible  for  the  peak  in  the  XPS  density  of  states  in  the 
vicinity  of  the  Fermi  energy.  The  highest  occupied  cluster  molecular  orbital 
8t1uUp)  indicated  in  Fig.  2  corresponds  to  Ti-bonding  combinations  of 
degenerate  and  A-j  p-band  Bloch  states  at  the  Fermi  energy,  and  the 

lowest  unoccupied  orbital  5 t^g  corresponds  to  o-nonbonding  p-band  Bloch 
states  just  above  the  Fermi  energy. 

As  will  be  shown  in  the  following  section  of  the  paper,  the  cluster 
t^uUp)  molecular  orbital  is  key  to  understanding  the  real -space  nature  of 
the  superconducting  state  of  aluminum.  Wave-function  contour  maps  of  this 
orbital,  plotted  in  the  (200) and  (110)  cyrstallographic  planes  up  to  second- 
nearest  neighbors,  are  shown  in  Figs.  4  and  5,  respectively.  The  solid  and 
dashed  contours  represent  positive  and  negative  regions,  respectively,  of 
the  molecular-orbital  wave  function.  Regions  of  net  overlap  between  contours 
of  the  same  sign  are  bonding,  whereas  regions  of  positive  and  negative  contours 
separated  by  nodes  are  anti  bonding.  The  most  striking  feature  of  Figs.  4  and  5 
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is  the  7r-bonding  character  of  the  A1  atomic  p  orbitals  along  the  "horizontal" 
direction  of  each  figure,  much  like  the  pi^ -bonding  character  of  C  atomic  p 
orbitals  around  the  ring  of  the  benzene  molecule  shown  in  Fig.  1.  Like  the 
planar  configuration  of  pn  atomic  orbitals  in  benzene,  the  pit-bonding  atomic 
orbitals  in  aluminum  form  two-dimensional  "layers"  as  shown  in  the  perspective 
contour  map  of  the  cluster  tjuUp)  molecular-orbital  wave  function  in  Fig.  6, 
which  is  the  three-dimensional  counterpart  of  Figs.  4  and  5.  However,  in 
contrast  to  benzene,  where  the  pit^  molecular  orbital  lies  well  below  the 
Fermi  energy,  i.e.,  it  is  the  "lowest"  occupied  orbital  in  Fig.  1,  the  pu-bonding 
molecular  orbital  of  aluminum  is  located  exactly  at  the  Fermi  energy,  i.e.,  it  is 
the  highest  occupied  cluster  molecular  orbital ,8t^u(ep) ,  in  Fig.  2.  Unlike 
benzene,  futhermore,  there  is  also  a  "one-dimensional"  a*-anti bonding  component 
between  the  p  orbitals  of  second-nearest-neighbor  A1  atoms  along  the  "vertical" 
direction  in  Figs.  4,5,  and  6,  which  effectively  isolates  the  pir-bonding 
atomic-orbital  components  between  nearest-neighbor  planes  of  atoms  and 
synergistical ly  promotes  covalent  bond  overlap  between  neighboring  atomic 
planes  to  form  the  layered  topology  of  the  molecular  orbital  shown  in  Fig.  6. 
Although  Figs.  4  through  6  imply  anisotropy  in  each  t1(j  Up)  molecular-orbital 
component,  the  triple  degeneracy  of  this  orbital  leads  to  equivalent  p^- 
bonding  and  po*-anti bonding  components  along  complementary  orthogonal  crystal 
directions  and  therefore  to  an  isotropic  three-dimensional  orbital  probability 
distribution.  As  the  cluster  size  is  systematically  increased  to  include  more 
shells  of  atoms  in  the  crystal,  each  pir-bonding  and  po*-anti bonding  molecular 
orbital  at  the  Fermi  energy  expands  spatially  and,  in  the  limit  of  a  very  large 
number  of  atoms,  forms  a  coherent  macroscopic  molecular-orbital  of  the  type 
shoym  schematically  in  the  (200)  crystal  plane  in  Fig.  7.  This  "real-space" 
character  of  aluminum  molecular  orbitals  at  the  Fermi  energy  cannot  be  extracted 
easily  from  the  conventional  "k-space"  Bloch  description  of  the  crystal  band 
structure,  except  in  hindsight  as  an  appropriate  combination  of  Zi ,  Zj,  and  kj 
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Bloch  eigenstates,  and  thus  has  eluded  solid-state  physics. 

III.  CLUSTER  ilOLECULAR  ORBITALS  AND  COOPER  PAIRS 

In  order  to  understand  how  cluster  molecular  orbitals  of  the  layered 
topology  shown  in  Fig.  6  can  be  precursors  to  Cooper  electron  pairing  and  the 
superconducting  state,  it  may  be  noted  first  that  each  layer  can  be  described 
by  a  bonding  orbital  wave-function  component,  y+  or  \p_,  where  the  plus  and 
minus  signs  designate  the  sign  of  the  wave  function.  Notice  also  in  Fig.  6 
that  along  the  direction  perpendicular  to  each  layer,  the  layer-orbital 
components  of  opposite  sign,  ip  and  y  ,  are  mutually  anti  bonding  and  cor¬ 
respond  approximately  to  a  damped  standing  sine  wave  having  nodes  that  coincide 
with  parallel  planes  of  atoms  and  a  wavelength  Xp  -  2d,  where  d  is 
the  interlayer  atomic  spacing.  In  other  words,  along  the  direction  perpendicular 
to  the  layer  orbitals,  the  wave  function  is  very  much  like  that  of  a  free 
electron  moving  back  and  forth  with  a  de  Broglie  momentum 

PF  «  h/\p  *  h/2d  (2) 

and  kinetic  energy  (the  Fermi  energy) 

eF  =  pF2/2m  *  h2/8md2.  (3) 

On  the  other  hand,  within  each  layer-orbital  component  of  Fig.  6,  the  wave 
function,  y+  or  y_,  corresponds  to  a  nodeless  standing  wave  of  wavelength 
X^  =  2L,  where  L  is  the  layer-orbital  diameter,  like  the  lowest-energy 
eigenstate  of  a  free  electron  confined  to  a  two-dimensional  constant  potential 
well.  With  increasing  cluster  size,  i.e.,  increasing  cluster  diameter  L,  the 
two-dimensional  component  of  de  Broglie  electron  momentum  parallel  to  the 
layers,  namely, 

PL  ■  h/XL  *  h/2L,  (4) 

vanishes 

lim  pL  *  0,  (5) 


L  -►  00 


-13- 


leaving  only  the  sine-wave-like  electron  motion  and  momentum  component  (3) 

perpendicular  to  the  layers  {see  Fig.  6).  However,  as  emphasized  in  Section 

II,  this  does  not  imply  a  preferred  direction  for  the  latter  component  of 

electron  momentum,  since  there  are  three  degenerate  molecular  orbitals  of 

the  topology  shown  in  Fig.  6  corresponding  to  three  orthogonal  crystal  directions 

,  % 

which  are  equally  accessible  at  the  Fermi  energy,  thereby  resulting  in  an 
effectively  isotropic  three-dimensional  sine-wave-like  electron  motion. 

From  one-electron  layer  orbitals,  *+  and  *  ,  of  the  type  shown  in  Fig.  6, 
it  is  possible  to  construct  an  electron-pair  wave  function  that  has  the  customary 
spin-paired  singlet  symmetry  of  a  Cooper  pair.  [1,2]  This  pair  wave  function 
can  be  written  as  the  product  of  a  symmetric  two-electron  space  wave  function 
*0,2)  =  *+(l)*_(2)  +  *_(l)*+(2)  (6) 

and  an  antisymmetric  two-electron  spin  wave  function 

S(l,2)  -  a(l)S(2)  -  3(l)ct(2),  (7) 

where  :s  and  i  are  functions  representing  "spin-up"  and  "spin-down"  electrons, 
respectively.  It  is  also  possible  to  construct  a  triplet  pair  wave  function 
from  the  ,+  and  orbitals,  although  we  shall  not  consider  this  configuration 
here.  Note  in  the  singlet  configuration  that  if  the  layer  molecular-orbital 
components  ;+  correspond  to  a  "spin-up"  (t)  electron,  then  the  layer  molecular- 
components  correspond  to  a  "spin-down"  (f)  electron  and,  vice  versa,  suggesting 
a  type  of  "anti ferromagnetic"  spin  polarization  of  the  alternating  layer  orbitals, 
i.e.,  a  "spin-density  wave"  (see  Fig.  6).  However,  the  net  spin  polarization 
per  atom  is  zero,  and  therefore  this  dt localized  "conduction-electron 
anti  ferromagnetism"  should  be  distinguished  from  conventional  localized 
anti  ferromagnetism  (see  further  discussions  in  Section  VI).  The  above  electron  pair 
wave  function  is  similar  in  form  to  the  ground-state  wavefunction  of  an  isolated 


hydrogen  (Hg)  molecule  in  the  valence-bond  description,  [55]  namely, 


,0,2)  *  ^a(l)*b(2)  ♦  *b(l)*a(2), 


(8) 
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where  <|>  and  are  the  atomic  (Is)  orbitals  of  the  dissociated  H  atoms. 

Q  D 

The  difference 

A  =  2els  -  E(d)  (9) 

between  the  dissociated-atom  energy  and  the  energy  E(d)  of  the  valence-bond 
state  (3)  for  the  equilibrium  internuclear  distance  d  defines  the  chemical  binding 
energy  of  the  H2  molecule.  A  stable  chemical  bond  (i.e.,  positive  value  of  A) 
will  not  occur  for  unless  the  Coulomb  and  exchange  repulsion  energies, 
respectively, 

URep  =  e2//|ipa(l )  |2i^b(2)  |2r12’1dvidv2,  (10) 

JRep  =  e2^V(1)V(2)^a^K(1)r12“ldvldv2  (11) 

of  the  two  electrons,  together  with  the  electrostatic  repulsion  of  the  two 
nuclei  (protons)  are  more  than  balanced  out  at  some  internuclear  distance  d 
by  the  Coulomb  and  exchange  attractions,  respectively, 

“attr  '  -«2--i*aO>i2’-,b',dvl  -«2/!»b(2)!2r2a-'dv2  (12) 

Jattr  ■  -d2S/ta*0)fb0)r,|)'’dv1  -e2S-‘ ^*(2) v&(2) dv2  (13) 

of  the  electrons  to  the  opposite  nuclei,  where 

s  =  V(1)V1)dvl  =  'V(2)V2)dv2  (14) 

is  the  atomic  overlap  integral. 

Analogously,  one  might  expect  the  energy  E(d)  for  the  "valence-bond-like" 
configuration  (6)  of  a  Cooper  pair  to  be  lower  than  the  energy  2t^  of  a 
"dissociated"  electron  pair  at  the  Fermi  energy,  leading  to  “quasimolecular" 
electron  pair  binding  of  the  type  envisioned  in  a  superconductor  by  Schafroth 
et  al_.  [56]  with  binding  energy 

A  =  2eF  -  E(d),  (15) 

if  the  Coulomb  and  exchange  repulsion  energies  of  the  electrons  are  more  than 
compensated  for  by  a  net  attraction  of  the  electrons  via  the  positive  ions 
of  the  lattice.  The  Coulomb  and  exchange  repulsions,  respectively, 

URep(cF)  “  e2//M1)i2|'*-(2)i2r12'1dv1dv2 


(16) 
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JRep(eF)  *  e2  I  ^+*0 )iP„*(2V+(2)ip_(l )r12'1dv1dv2  (17) 

of  a  pair  of  electrons  correlated  at  the  Fermi  energy  over  degenerate, 
mutually  antibonding  layer  orbitals  of  the  type  shown  in  Fig.  6  are  already 
largely  screened  out  over  distances  greater  than  the  interlayer  spacing  d  by 
the  positive  ions  located  in  the  nodal  planes  between  the  layer  orbitals. 


Then  the  quantities  (16)  and  (17)  reduce  approximately  to 

URep(eF)  ^  (eZ/d)/|i|,+(1)|2dv1/|ip_(2)|2dv2  =  e2/d  (18) 

JRep(eF)  ^  (e2/d)/^*n)^0)dv1/w.*(2),i,+  (2)dv2  =  0  (19) 

since  the  layer  orbitals  are  normalized,  i.e., 

/k+(1)|2dv1/|^.(2)|2dv2  =  1  (20) 

and  mutually  antibonding,  i.e.,  their  overlap  integral 

S+  =  Sant1_UF)  »  .%(l)v_(l)dv1  =  ;*+(2)*_(2)dv2  •  0.  (21) 

~  bond 


Therefore,  the  screened  electrostatic  repulsion  of  two  such  electrons  at  a 


relative  distance  r  can  be  described  approximately  by  the  effective  potential 


VC<r)  v  u«*pUF> 


e  /d  for  r  <  d 
0  for  r  >  d 


where  d 


2A  for  aluminum.  The  Coulomb  attraction 
UAttr(sF}  =  -e2/«*4i+(1)  |2r1“1dv1  -e2/|^_(2)  |2r2‘] 


dv. 


(22) 

(23) 


of  each  layer  orbital  to  the  intervening  positive  ions  is  already  accounted  for 
in  the  potential  which  determines  the  molecular-orbital  energy  2£p  of  the 
"dissociated"  electron  pair,  while  the  exchange  attraction 

JAttr(£F)  =  -e2S+.^+(l)’|/.(l)r1-1dv1  -e2S+/^(2)v_(2)r2'1dv2  (24) 

vanishes  because  the  overlap  integral  (21)  between  the  mutually  anti  bonding 
layer  orbitals  p+  and  t{/_  is  zero. 

Because  of  the  molecular-orbital  degeneracy  at  the  Fermi  energy,  the 
electron  pair  described  by  the  layer  orbitals  -p+  and  p_  is  subject  to  the 
Jahn-Teller  effect.  [57]  Therefore,  the  electron  pair  can,  in  principle,  lower 
its  energy  from  the  value  2tp  by  the  amount  (15)  by  coupling  with  the  lattice 
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through  displacements  of  the  lattice  ions.  The  "driving  force"  for  the  Jahn- 
Teller  coupling  is  the  anisotropy  of  the  degenerate  :p+  and  -j/_  layer  orbitals, 
which  tends  to  attract  the  nearest-neighbor  lattice  ions  along  the  pTr-bonding 
directions  parallel  to  the  layers  and  to  repel  second-nearest-neighbor  ions 
along  the  pc*-anti bonding  direction  perpendicular  to  the  layers  (see  Figs.  4 
through  7).  An  instantaneous  displacement's  of  a  lattice  ion  with  respect 
to  the  neighboring  layer  orbitals  ’J/+(1)  and  p_(2)  induces  an  attractive 
electrostatic  dipole  potential  of  the  form 

UAttr)(eF)  =  -e2'b+0)|2r1“26dv1  -e2/jj,_(2)  |2r2  “25dv2.  (25) 

For  the  range  of  ion  displacements 

6  \  (m/M)3d  for  1/4  <  B  <  1/2  (26) 

1  /4 

between  the  amplitude  6q  ^  (m/M)  A  of  the  zero-point  oscillation  of  the 

1  /2 

lattice  and  the  amplitude  5c  w  (m/M)  '  d  of  a  classical  oscillator  (where 
m  is  the  electron  mass  and  M  is  the  ion  mass),  it  is  unfavorable  for  the 
electron-lattice  system  to  remove  the  orbital  degeneracy  at  the  Fermi  energy 
and  lower  the  electronic  energy  through  a  static  lattice  distortion.  [57] 
Instead,  the  -p+  and  p_  orbitals  and  lattice  ion  displacements  o  are  coupled 
dynamically,  and  the  range  of  3  values  in  (26)  corresponds  to  the  so-called 
dynamic  Jahn-Teller  effect.  [57]  From  the  SCF-Xa-SW  cluster  molecular-orbital 
contour  map  in  Fig.  6,  which  is  a  direct  measure  of  the  average  spacing  d  of 

O 

the  "compressed"  layer  orbital  components  and  a  value  of  d  v  1 .3A 
appropriate  for  bulk  crystalline  aluminum  is  obtained,  giving  limiting  values 

O  0 

of  ^  0.1  A  and  %  0.01  A  for  aluminum. 

0  c 

The  Jahn-Teller  theorem  [57]  predicts  only  the  occurrence  of  a  coupling 
between  orbitally  degenerate  electrons  at  the  Fermi  energy  and  lattice-ion 
displacements  and  not  their  actual  magnitude.  Nevertheless,  it  is  reasonable  to 
expect  that  the  Jahn-Teller-induced  lattice-ion  displacements  5  are  proportional 
to  the  overlap  integral 
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Sbond(eF)  =  *  /*_a(l)H>_b(l)dv1  (27) 

between  nearest-neighbor  atomic-orbital  components  a  and  b  of  the  molecular- 
orbitals  ip+  and  ip_  at  the  Fermi  energy  along  the  bonding  direction,  e.g.  the 
pir-bond  overlap  within  the  ^-bonding  layers  of  aluminum  shown  in  Figs.  4  through  7. 
Then  one  can  write  approximately 

S  '  W£F)d-  <“> 

where  d  is  of  the  order  of  the  interatomic  distance.  The  bond  overlap  integral, 
Sbond^F^’  sh0u^  1)6  distinguished  from  the  zero  overlap  integral  $  .  Up) 

defined  in  (21)  between  mutually  anti  bonding  layer  molecular-orbita*?on<^ 
components  in  Fig.  6.  In  the  limit  of  relatively  small  Sjj0n<jUp)*  it  is  also 
reasonable  to  expect  that  the  lattice-ion  displacements  6  should  approach  the 
classical  limit  cc  ^  (m/M)  '  d  mentioned  above,  whereas  in  the  limit  of 
reasonably  large  ^bond^F^’  5  should  be  bounded  by  the  quantum-mechanical  limit 
ig  (m/M)1//4d.  In  other  words,  comparing  expressions  (28)  and(26),  there  is 
an  implied  relation 


(m/M)~  ^  Sbond(“;F)  (29) 

between  the  Jahn-Teller  coupling  parameter  £  and  the  bond  overlap  at  the 

Fermi  energy.  For  aluminum,  where  the  pr-bond  overlap  at  the  Fermi-energy,  as 
evaluated  from  the  molecular-orbital  wave  functions  mapped  in  Figs.  4  through  6, 
is  small,  expression  (29)  implies  a  value  of  3  ^  1/2.  This  valje  of  3  defines 
the  limit  of  "weak"  Jahn-Teller  coupling,  in  which  Jahn-Tel ler-induced  lattice 
ion  displacements  5  defined  in  (26)  are  inversely  proportional  to  the  square 
root  of  the  ion  mass,  and  corresponds  to  the  normal  "isotope  effect"  in  super¬ 
conductivity.  Indeed,  as  will  be  elaborated  on  in  this  paper,  the  dynamic  Jahn- 
Teller  coupling  of  electrons  and  lattice  is  the  real-space  molecular-orbital 
analogue  of  the  coupling  of  electrons  and  virtual  phonons  responsible  for  Cooper 
pair  formation  in  the  BCS  theory  of  superconductivity.  [1,2]  However,  in  the  present 
model,  electron-lattice  coupling.  Cooper  pair  formation,  and  superconductivity 
will  not  occur  unless  there  is  a  prescribed  molecular-orbital  topology  at  the 
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Fermi  energy,  such  as  that  shown  for  aluminum  in  Fig.  6. 

In  the  limit  of  dynamic  Jahn-Teller  coupling,  all  components  of  the 
degenerate  molecular  orbital  at  the  Fermi  energy,  including  the  one  shown  in 
Fig.  6  and  equivalent  ones  pointing  along  orthogonal  crystal  directions, 
are  dynamically  accessible,  resulting  in  an  effectively  isotropic  attractive 
potential  (25)  of  two  electrons.  From  (25)  and  (26),  this  interaction  can  be 


written  as  a  simple  pair  potential 

VA(r)  %  UA2t  (eF}  ^{'e25/d2  ^  -e2(m/M)S/d^  -e2/£  for  r  <  l 


0  for  r  >  l 


where  r  is  the  relative  distance  of  the  two  electrons  and 

,3 


(30) 

(31) 


l  \  (M/m)  d 

is  the  effective  distance  and  ultimate  cluster  radius  over  which  the  dynamic 

Jahn-Teller  effect  and  attractive  pair  potential  (30)  are  operative.  In 

1  ° 

aluminum,  where  3  y,  l  ^  400A.  Combining  (22)  and  (30),  one  finds  the  total 
effective  potential  between  the  valence-bond-like  pair  of  electrons  at  the 
cermi  energy  to  be 


V(r) 


Vc(r)  ♦  VA(r) 


'{ 


e2/d  -e|/ i  for  r  <  d 
-e  /i  for  d  <  r 
0  for  r  >  l 


(32) 


0  3 

Thus,  over  a  distance  d  <  r  <  l  spanning  many  atoms  (2A  <  r  <  400A  in  aluminum), 
the  effective  potential  of  two  electrons  at  the  Fermi  energy,  induced  by  the 
dynamic  Jahn-Teller  effect,  is  weakly  attractive  relative  to  the  average 
equilibrium  lattice  potential.  As  will  be  argued  below,  this  attractive  potential 
leads  to  Cooper-Schafroth  [2,56]  "quasimolecular"  electron-pair  binding  of  the 
type  prerequisite  for  superconducti vity . 

While  a  value  of  3  A-  1/2  in  (26),  (29),  (30),  and  (31)  is  appropriate  for 
aluminum  and  other  "weakly"  Jahn-Teller  coupled  p^r-electron  superconductors 
exhibiting  the  normal  isotope  effect, values  of  3  <  1/2  are  more  appropriate 
for  transition-metal  superconductors.  As  will  be  discussed  in  Section  IV, 
transition-metal  superconductors  are  characterized  by  degenerate  molecular 
orbitals  at  the  Fermi  energy  of  "tubular"  topology  (see  Figs.  11  through  14)  arising 
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from  d5-bond  overlap  and  of  "layered"  topology  (see  Fig.  16)  arising  from  d^2 
planar  overlap.  These  orbitals  are  argued  to  be  responsible  for  Jahn-Teller 
coupling  and  Cooper-Schafroth  pairing  in  transition-metal  superconductors, just  as 
pr  orbitals  are  responsible  for  these  effects  in  aluminum.  In  general ,  however, 
d6-bond  overlap  in  transition  metals  is  greater  than  pTr-bond  overlap  in  non¬ 
transition  metals  such  as  >  luminum,  implying  from  (26),  (28),  and  (29)  stronger 
Jahn-Teller  coupling  and  values  of  3  <  1/2  for  transition  metal  superconductors 
consistent  with  deviations  from  the  ideal  isotope  effect  (see  Section  IV). 

If  bond  overlap  at  the  Fermi  energy  is  large  enough  to  yield  a  value  of 
3  <  1/4,  the  Jahn-Teller  coupling  can  result  in  lattice  ion  displacements  6 
greater  than  the  zero-point  amplitude  5g  ^  (m/M)^4d,  i.e., 

o  ^  (m/M)Sd  for  3  <  1/4.  (33) 

Then  the  electrons  at  the  Fermi  energy  can  lower  their  energy  through  a 
periodic  three-dimensional  static  Jahn-Teller  lattice  distortion  [57]  (or 
alternatively  a  "one-dimensional"  Peierls  distortion  [58])  which  removes  or 
lowers  the  orbital  degeneracy  and  opens  up  an  energy  gap  around  the  Fermi 
energy.  This  process  can  compete  with  the  electronic  energy  reduction  effected 
by  the  dynamic  Jahn-Teller  effect,  Cooper-Schafroth  pairing,  and  the  opening  up 
of  the  superconducting  state  energy  gap  A  defined  in  Eq.  (15) .  Indeed  structural 
phase  transitions  corresponding  to  static  Jahn-Teller  distortions  of  the  crystal 
lattice  from  cubic  to  lower  (e.g.,  tetragonal)  symnetry  are  well  known  to 
occur  in  certain  high-temperature  transi tion-metal -compound  superconductors 
(e.g.,  A15  compounds),  [58]  and  such  "lattice  instabilities"  are  suspected  to 
reduce  the  superconducting-state  transition  temperature  Tc  otherwise  attain¬ 
able.  [58,3]  In  any  case,  the  conjecture  that  the  interactions  responsible  for 
these  structural  transformations  and  high-temperature  superconductivity  are 
similar  [58]  is  supported  by  the  present  theoretical  model. 

Finally,  when  the  degeneracy  of  the  spatially  extended  bonding  molecular 
orbitals  at  the  Fermi  energy  is  removed  or  lowered  by  a  static  Jahn-Teller 
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distortion,  Peierls  distortion,  or  intrinsically  low  (non-cubic)  crystal 
symmetry,  it  is  still  possible,  within  the  framework  of  this  model,  for 
Cooper-Schafroth  electron  pairing  and  superconductivity  to  occur  via 
dynamic  "pseudo-Jahn-Teller"  electron-lattice  coupling  [57]  triggered  by 
dynamical  fluctuations  in  electron  occupancy  or  "mixing"  of  the  nearly 
degenerate  molecular  orbitals.  8ecause  of  the  different  spatial  orientations 
of  these  orbitals,  this  process  can  involve  fluctuations  in  orbital  charge 
distribution  around  the  Fermi  energy  and  therefore  fluctuations  5  in  the 
relative  instantaneous  locations  of  bond  charge  and  lattice-ion  positions 
at  the  Fermi  energy  approaching  the  interorbital  spacing  d.  Such  large 
relative  displacements  of  electrons  and  nuclei  imply,  according  to  (33), 
a  value  of  3  approaching  zero.  Since  fluctuations  in  electron  occupancy  or 
electronic  "excitations"  between  nearly  degenerate  molecular  orbitals  around 
the  Fermi  energy  are  associated  with  the  dynamic  "quasi -Jahn-Tel ler"  effect,  [57] 
this  process  can  be  viewed  as  having  a  purely  electronic  (i.e.,  "non-phonon") 
component.  However,  this  "excitational "  mechanism  is  dependent  on  the  same 
molecular-orbital  topologies  prerequisite  to  the  above-described  conventional 
dynamic  Jahn-Tel ler  (electron-phonon)  coupling  mechanism  and  thus  should  not 
be  confused  with  "exci tonic"  mechanisms  for  superconductivity  proposed  by  other 
workers.  [10,11]  Like  the  ordinary  dynamic  Jahn-Teller  effect,  the  dynamic 
"quasi-Jahn-Teller"  effect  can  produce  an  attractive  electron  pair  potential 
of  the  form  (30),  although  since  5  n,  d  and  3^0,  the  potential  is  of  short- 
range  character  and  is  largely  cancelled  out  by  the  screened  repulsive  Coulomb 
potential  (22).  It  is  possible  that  the  latter  mechanism  is  primarily  responsible 
for  Cooper-Schafroth  pairing  and  superconductivity  in  certain  non-cubic  transition 
models  (e.g.,  Zr  and  Ru),  transition-metal  compounds  (e.g.,  A15  and  "layered" 
compounds),  and  amorphous  transition-metal  alloys  (e.g.,  CuZr),  which  exhibit 
vanishing  (3  ^  0)  or  anomalous  isotope  effects  and  where  the  molecular  orbitals 
at  the  Fermi  energy  responsible  for  superconductivity  are  dc-bonding  "tubular" 
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orbitals  of  the  types  shown  in  Figs.  11  through  14  (see  Section  IV). 

In  order  to  make  a  more  formal  connection  between  this  model  and  the 
BCS  theory  of  superconductivity,  [1]  one  can  make  use  of  simple  physical 
arguments  recently  proposed  by  Weisskopf  [59]  for  the  formation  of  Cooper  pairs 
in  the  BCS  theory.  In  Weisskopf's  model  for  a  superconductor,  an  electron 
moving  at  the  Fermi  energy  with  momentum  (2)  and  speed 

vg  *  pF/m  3  h/2md  (34) 

transfers  some  of  its  momentum  to  neighboring  ions  of  the  lattice,  causing 
ion  displacements  6  toward  the  electron  path  of  magnitude  (26)  in  the  limit 
3  'v  -a—.  However,  the  ions  do  not  remain  in  their  displaced  positions  in¬ 
definitely  but  return  to  their  equilibrium  lattice  positions  in  a  time  of 
order  ,  where 

a»0  h/2(mM)1/2d2  (35) 

is  the  approximate  expression  for  the  lattice  Debye  frequency  derived  by 
Weisskopf.  [59]  Thus,  an  ion  displacement  toward  the  electron  path  can  be 
viewed  as  extending  over  a  distance 

1  ^  ve^u,D  (36) 

behind  the  electron  and  producing  over  this  distance  an  attractive  potential 

for  a  second  electron  moving  "head-on"  in  the  opposite  direction.  This  is 

Weisskopf's  simple  physical  picture  of  the  emission  of  a  virtual  phonon  by 

one  electron  and  the  absorption  of  the  phonon  by  a  second  electron  responsible 

for  electron-electron  attraction  and  Cooper  pair  formation  in  the  BCS  theory.  [1,2] 

Substitution  of  (34)  and  (35)  into  (36)  yields  the  expression  (31)  for  l 

derived  above  in  the  limit  3  ^  ^  from  ttie  dynamic  Jahn-Teller  effect,  so  that  the 

effective  attractive  electron-pair  potential  can  be  written  in  the  form  (30), 

which  Weisskopf  views  as  being  contained  in  a  "long  thin  tube"  of  length  i 

and  diameter  d  while  acting  only  on  electrons  of  opposite  linear  momentum, 

opposite  spin,  and  in  a  relative  "S-state"  (relative  angular  momentum  L  *  0). 

Thus,  Weisskopf's  model  and  the  present  molecular-orbital  model  are  in 
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basic  agreement  as  to  the  general  form  of  the  attractive  potential  responsible 
for  Cooper  pairing  in  a  superconductor  in  the  limit  3  %  and  both  models 
lead  to  simple  "real -space"  pictures  of  the  electron-lattice  coupling  which 
are  analogous  to  the  electron-phonon  coupling  mechanism  in  the  BCS  theory. 
However,  Weisskopf’s  intuitive  picture  of  the  electron-lattice  coupling  in 
a  superconductor,  like  the  original  BCS  theory,  does  not  address  the  question 
of  which  materials  are  superconductors  and  which  are  not,  and  it  does  not 
account  for  superconductors  with  vanishing  isotope  effect  (3^0).  On  the 
other  hand,  the  molecular-orbital  model  suggests  that  superconductivity  will 
not  occur  unless  there  are  certain  types  of  spatially  extended  orbital  topology 
(e.g.,  like  that  in  Fig.  6)  at  the  Fermi  energy,which  promote  Cooper-Schafroth 
electron  pairing  via  the  dynamic  Jahn-Teller  effect  for  3  \  j  and  via  the 
dynamic  "quasi -Jahn-Tel ler”  effect  for  3  \.  0. 

To  determine  the  quasi  molecular  binding  energy  (15)  of  a  Cooper-Schafroth 
electron  pair  and  the  superconducting  transition  temperature  T_  from  the  present 
theoretical  model,  one  can  follow  Weisskopf's  [59]  procedure  and  solve 
Scnrodinger 's  equation 

[p2/2u  +  V(r)]v(r)  =  E^(r)  (37) 

for  two  electrons  of  relative  distance  r,  relative  linear  momentum  p,  total 
linear  momentum  P  =  0,  relative  angular  momentum  L  =  0,  and  reduced  mass 
u  =  m/2,  subject  to  the  potential  (32).  Because  this  potential  is  attractive 

o  o 

over  a  distance  d  <  r  <  l  (2A  <  r  <  400A  in  A1 )  spanning  many  atoms  and 

O 

vanishes  for  r  >  i,  the  repulsive  component  for  0  <  r  <  d  (0  <  r  <  2A  in  A1 ) 
should  not  strongly  influence  the  binding  energy  and  large  spatial  extent  of 
the  pair  wave  function.  Thus  one  can  solve  Eq.  (37)  to  good  approximation  only 
for  the  attractive  component  (30)  of  the  potential,  subject  to  the  boundary 
condition  v( d)  =  0. 

The  Schrodinger  equation  (37)  is  effectively  a  one-dimensional  wave 
equation  whose  solution  p(r)  is  actually  r  times  the  three-dimensional  valence- 
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bond-like  wave  function  (6).  Since  the  latter  wave  function  is  expressed  in 
terms  of  the  cluster  layer  orbitals  tp+  and  -,j/_  shown  in  Fig.  6,  which  have  damped 
standing-sine-wave-like  character  along  r,  it  is  appropriate  to  expand  the 
solution  of  (37),  namely 

'+>( r )  =  ,  £  a ( p '  )tp( p '  ;r) ,  (38) 

P  >  PF 

in  the  set  of  sinusoidal  basis  functions 

’•Hp;**)  s  ( 2/R) 1  /2si n  (pr/h)  (39) 

normalized  to  a  large  sphere  of  radius  R  >  i,  as  given  by  (31).  This  sphere 

radius  defines  the  ultimate  cluster  size  necessary  to  describe  the  full  spatial 

0 

extent  of  the  quasimolecular  electron-pair  wave  function  (e.g.,  R  >  400A  in  Al). 

The  summation  in  (38)  excludes  the  momenta  p  <  Pp  corresponding  to  other 
electrons  Delow  the  Fermi  energy,  and  for  a  pair  wave  function  of  large  spatial 
extent  the  coefficients  a(p)  will  be  large  only  for  p  ^  Pp  ^  h/2d  [see  Eq.  (2)]. 
Note  that  for  such  values  of  momentum  the  basis  functions  (39)  satisfy  the 
boundary  condition  j(p;d)  =  0.  It  should  be  emphasized  that,  while  the  basis 
functions  (39)  are  identical  to  the  free-electron  basis  set  adopted  by  Weisskopf 
[59 J  ana  a-e  equivalent  ot  the  plane-wave  basis  of  the  BCS  theory,  [1]  their  simple 
sinusoidal  form  is  actually  a  by-product  of  the  layered  topology  of  the  composite 
bonding-antibonding  molecular  orbitals  at  the  Fermi  energy  shown  in  Fig.  6.  Thus 
the  following  derivation  of  the  Cooper  pair  binding  energy  and  transition 
temperature  for  superconductivity  using  this  basis  set  is  strictly  valid 
only  for  substances  that  have  the  prerequisite  molecular-orbital  topology  at 
the  Fermi  energy. 

Substitution  of  (38)  into  (37),  multiplication  of  both  sides  of  the 
resulting  equation  by  (39),  and  integration  over  the  sphere  of  radius  R  yields 
a(p)(p2/m  -  E)  =  -I  a(p’ )/R4p;r)V(r)v(p' ;r)dr.  (40) 

p '  >  pp  0 
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Using  the  identity 

v(p;r)4>(p' ;r)  *  (l/R){cos[(p-p')r/fi]  -  cos[ (p  +  p')r/R]},  (41) 

and  the  fact  that  the  largest  terms  in  the  summation  of  (40)  occur  for  momenta 
p'  x  p  -v,  pF  ^  h/2d,  for  the  potential  given  in  (30),  the  integral  in  (40) 
reduces  approximately  to 

rRib(o-r)\J(r'm(o' -rldr  J'd/2')  for  P  “  P*  <  F i/4 

'n  4'(P’r,V(r/1"(P  ,r)dr  0  for  p  -  p1  >  fi/Jl  (42) 


or  substituting  (31)  in  the  right  side  of  (42),  to 


rK 


j  ’^(p*»r ) V(r)u>(p 1  ;r)dr  'u 
0 


t 


-(e  /R)[l 
0 


■(m/M)S]  for  p  -  p 1 
for  p  -  p 1 


<  R/(M/m)?d 
>  R/(M/m)Bd, 


(43) 


which  is  an  attractive  pair  pseudopotential  for  l  >  d  or  g  >  0.  Since  the  electron 
pair  binding  energy  is  to  be  determined  with  respect  ot  the  Fermi  energy  of 
otherwise  independent  electrons,  it  is  convenient  to  introduce  the  new  energy 
variable 


e  =  (p2/m)  -  2cp  =  (p2  -  Pp2)/m 

and  to  rewrite  Eq.  (40)  in  terms  of  s  and  the  quasimolecular  pair  binding 
energy  A  defined  in  Eq.  (15).  Thus,  in  place  of  (40),  one  can  write 
a(e)(e  +  A)  =  -£  a(e‘  )/R<p(e  ;r)  V(r)^(s  *  ;r)dr, 

^  1  ^  n  0 


(44) 


(45) 


where  the  summation  excludes  the  energies  corresponding  to  electron  states 
below  the  Fermi  energy,  since  such  states  do  not  contribute  to  Cooper  pairing. 
The  energy  interval  t  -  e1  that  corresponds  to  the  maximum  momentum  difference 
P  -  p1  v  h/i  in  (42)  is,  using  expressions  (2)  and  (31), 

ie  ■  (e  •  £,)™x  '  (2pF/")<p  -  p,)n«x  (46) 

'u  (h/md)(fiM) 

■v  h2/2-rm(M/m)3d2. 

Thus  the  integral  in  (45)  reduces  to  a  pseudopotential  expression  analogous 
to  (43),  namely. 
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%  t  1 1'.  I  * 

Substitution  of  (47)  in  (45)  and  conversion  of  the  summation  in  the  latter 

equation  to  an  integration  yields 

a(e)(e  +  A)  ^  N(0)(eVR)[1  -  (m/M)a]/a(e *  )de ' , 

0 


(47) 


(48) 

(49) 


where  ;r 

N(0)  %  (R/2ir) (m/hpp)  %  2mRd/h 
is  the  number  of  sinusoidal  basis  functions  (39)  of  de  Broglie  wavelength 
\p  =  2d  [see  Eq.  (2)]  per  molecular  orbital  (e.g.,per  layered  molecular 
orbital  in  Fig.  6)  per  unit  energy  dc  around  the  Fermi  energy  gp,  normalized 
to  the  sphere  of  radius  R.  Expression  (49)  defines  the  molecular-orbital 

"density  of  states"  at  the  Fermi  energy.  Substitution  of  (49)  into  (48)  yields 

Ac  + 


a(s)(e  +  a)  %  (2me2d/h2)[l  -  (m/M)®].'a(e')de' 

0 

a-  (2me  d/n  )[1  -  (m/M)^]/a (e 1  )de‘ 


(50) 


in  which  the  R-dependence  of  the  orbital  density  of  states  (49)  and  the  R- 

dependence  of  the  pseudopotential  (e  /R)[l  -  (m/M)3]  are  cancelled  out  through 

their  product  in  (48).  Furthermore,  because  the  coefficients  a(p')  in  (40) 

are  significantly  different  from  zero  only  for  p'  close  to  pF,  the  coefficients 

a(c’)  in  (SO)  are  significant  only  for  c’  «  Ae,  thus  allowing  one  to  use  Ae 

as  the  upper  limit  of  integration  in  (50).  Since  the  resulting  integral  in 

(50)  is  independent  of  s,  it  can  be  set  equal  to  the  constant 

Ae 

C  =  /a ( e ' ) dc 1 ,  (51) 

0 

and  Eq.  (50)  can  be  solved  for  a(s),  giving 

a(e)  x  (2me2d/h2)[l  -  (m/M)3]C(e  +  A)'1.  (52) 

Substituting  (52)  for  e  =  e1  in  (51)  and  cancelling  out  the  constant  C,  one 
finds 
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1  %  (2me2d/h2)[l  -  (m/M)e]/(e'  +  A)_1de' 

0 

^  (2me2d/h2)[l  -  (m/M)S]£n[(Ae  +  A)/A]  (53) 

^  (2me2d/h2)[1  -  (m/M)6]£n(Ae/A) , 

where  the  final  expression  (53)  results  from  the  fact  that  the  quasimolecular 
electron-pair  binding  energy  A  should  be  significantly  less  than  the  maximum 
energy  interval  Ae  defined  in  (46).  Substituting  (46)  in  (53)  and  solving 
for  A,  one  obtains 

A  %  [h2(m/M)3/2™d2]exp{-h2/2me2d[l  -  (m/M)3]}.  (54) 

Using  the  established  relation  between  the  superconducting  energy  gap  and 
transition  temperature,  namely,  [1,60,61] 

A  'v  2kgT  (or,  more  precisely,  kgT  =  0.57  A),  (55) 

where  kg  is  Boltzmann's  constant,  expression  (54)  can  be  reformulated  in  terms 
of  T  ,  giving 

kcT  [h2(m/M)3/4-md2]exp[-h2/{2me2d[l  -  (m/M)2]}]  (56) 

In  turn,  expression  (56)  can  be  put  into  the  form 

kgTc  ^  Ruic  exp[-l/(X  -  y*)]  (57) 

where 

hu)c  =  h2(m/M)3/4-nnd2  (58) 

X  =  2me2d/h2  =  N(0)(e2/R)  =  N(0)V  (59) 

and 

u*  =  (2me2d/h2)(m/M)6  =  N(0) (e2/R) (m/M)3  =  N(0)U.  (60) 

The  molecular-orbital  "density  of  states"  at  the  Fermi  energy,  N(0)  defined  in 
Eq.  (49),  has  been  incorporated  into  expressions  (59)  and  (60). 

Expression  (56)  reformulated  as  (57)  is  the  product  of  a  preexponential 
function  of  the  "cut-off  frequency"  wc  defined  in  (58)  and  an  exponential 
function  of  the  attractive  electron  pair  interaction  parameter  \  (or  attractive 
electron  pair  pseudopotential  V)  defined  in  (59)  and  the  repulsive  screened 
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Coulomb  interaction  parameter  u*  (or  Coulomb  pseudopotential  U)  defined  in  (60). 
Thus  (56)  and  (57),  derived  directly  from  molecular-orbital  theory,  are  in 
the  form  of  a  generalized  BCS-like  expression  for  the  superconducting  transition 
temperature  which  allows  for  values  of  the  Jahn-Teller  coupling  parameter 
S  <  j  [see  (26)  and  (33)]  and  consequently  is  applicable  to  all  superconducting 
materials,  including  tra.  4tion  metals  (see  Section  V)  where  the  simple  M"^2 
dependence  of  Tc  predicted  by  the  original  BCS  theory  [1]  is  not  confirmed  by 
experiment  [60]. 

To  relate  the  present  molecular-orbital  formulation  to  the  original  BCS 

theory,  [1]  we  consider  the  reduction  of  expressions  (56)  and  (57)  for  the  value 
1  8 

3  'v  -.  In  this  limit,  the  factor  (m/M)  in  (56)  is  negligible  compared  with 
unity  and  u*  <<  X  in  (57).  Furthermore,  for  3  ^  the  cut-off  frequency  u>c, 
defined  in  (58),  reduces  to  the  lattice  Debye  frequency  defined  in  (35).  There¬ 
fore,  in  place  of  (56)  and  (57),  one  can  write 

lim  kaT  [h2/4Tr(mM)1^d2]exp(-h2/2me2d) ,  (61) 

3  ■+  1/2  °  c 


and,  using  (35)  and  (59), 

lim  kgT.  %  fiuig  exp(-l/X) 

3  ^2  -v  hu)Q  exp[-l/N(0)V] 

Finally,  since  the  Debye  frequency  u,. 


“D 


through  the  formula 
kB0D  =  ^D’ 


is  related  to  the  Debye  temperature 


(62) 

r*v 

C'D 

(63) 


the  last  expression  in  (62)  reduces  to 
lim  T.  n,  0  exp[-l/N(0)V], 

a  -  1/2  c  u 

which  is  identical  in  formto  the  original  BCS  expression  (1)  for  T  . 

The  original  BCS  formula  (1)  for  the  superconducting  transition  temperature 
has  .been  difficult  to  apply  from  first  principles  to  real  materials  because  it 
presumes  a  knowledge  of  the  band-structure  density  of  states  at  the  Fermi  energy, 
N(0),  and  the  attractive  electron-phonon  coupling  potential  V.  These  two  quantities 
are  often  not  known  to  high  accuracy,  especially  V,  and  therefore 


must  be  estimated  semiempirically  to  correlate  with  experimentally  determined 

values  of  T  .  Furthermore,  for  certain  transition  metals,  their  alloys  and 
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compounds,  the  simple  dependence  of  T  on  M  '  predicted  by  the  BCS  formula 

v» 

(1)  [see  also  expression  (61)]  is  not  confirmed  experimentally  by  the  "isotope 
effect. "[60]  In  attempts  to  explain  the  variation  in  isotope  effect  among 
superconductors ,  without  invoking  alternative  mechanisms  for  Cooper  pairing, 
the  original  BCS  theory  has  been  generalized  by  various  workers  [60,64]  to 
include  a  more  accurate  treatment  of  the  combined  effects  of  Coulomb  electron- 
electron  repulsion  and  phonon-induced  electron-electron  attraction  in  determining 
the  net  potential  V  in  expression  (1).  The  formulae  for  Tc  in  these  generalized 
versions  of  BCS  theory  are  in  the  form  of  expression  (57),  although  the  specific 
definitions  of  the  parameters  fituc,  <\,  and  u*  differ  from  expressions  (53) 
through  (60)  and  vary  somewhat  with  the  investigator.  [60,64]  Exciton 
models  [10,11]  and  other  "non-phonon"  mechansism  [65-67]  for  superconductivity 
also  assume  a  generalized  BCS-like  formula  (57)  for  T  ,  in  which  FLC  is  the 
maximum  energy  of  the  electronic  excitation,  \  is  the  electron-exci ton  coupling 
constant,  and  u*  is  the  screened  Coulomb  repulsion.  Thus  far,  there  is  no 
definitive  experimental  evidence  for  the  existence  of  superconductors  based  on 
excitonic  or  other  non-phonon  mechanisms  of  the  types  described  in  Refs. 

10,  11 ,  65,  66,  and  67. 

The  present  molecular-orbi tal  model  for  Cooper-Schafroth  pairing  and 
superconductivity  complements  BCS  theory  and  clarifies  the  issue  of  non¬ 
phonon  mechanisms  for  superconductivity  in  the  following  respects: 

1.  It  permits  one  to  assess  which  materials  are  likely  to  be  superconductors 
and  which  are  not  on  the  basis  of  simple  criteria  for  extended  molecular- 
orbital  topology  at  the  Fermi  energy,  as  deduced  from  SCF-Xci-SW  cluster 
molecular-orbital  calculations. 

It  permits  the  approximate  calculation  of  the  superconducting  transition 
temperature  via  Eq.  (56)  directly  in  terms  of  two  simple  parameters  derived 


2. 
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from  first  principles  SCF-Xa-SW  cluster  molecular-orbital  studies,  namely, 
a  local  distance  parameter  d  of  the  order  of  magnitude  of  the  interatomic 
distance  and  a  Jahn-Teller  coupling  parameter  8,  instead  of  in  terms  of 
the  electronic  density  of  states  N{0),  electron-phonon  pseudopotential  V, 
and  Coulomb  pseudopotential  U. 

,  * 

3.  While  the  present  molecular-orbital  model  reduces  to  a  BCS-like  theory 

for  weakly  Jahn-Teller-coupled  pir-bonded  superconductors  such  as  aluminum, 

where  S  ^  it  is  also  readily  applicable  to  more  strongly  Jahn-Teller- 
2 

coupled  d6-bonded  transition-metal  superconductors,  including  compounds 

and  alloys,  that  exhibit  an  anomalous  "isotope  effect,"  i.e.,  where 

8  <  -  (see  Section  IV).  In  this  respect,  Eq.  (56)  written  in  the  form  (57) 

2 

already  includes,  without  further  generalizations  of  the  types  described 
in  Refs.  60-64,  the  effects  of  screened  Coulomb  repulsion  through  the 
dependence  of  both  the  quantity  u*  [defined  in  (60)]  and  pre-exponential 
energy  h^c  [defined  in  (58)]  on  the  factor  (m/M)  ,  which  increases  rapidly 
with  decreasing  2. 

4.  For  1/4  <  B  <  1/2,  where  the  relative  displacement  (26)  of  the  electron 

pair  and  lattice  ions  is  less  than  the  zero-point  vibrational  amplitude 
1  /4 

5q  (m/M)  d  of  the  lattice,  dynamic  Jahn-Teller  coupling  of  the  electrons 
and  lattice  ions,  driven  by  the  degeneracy  of  the  spatially  extended 
(e.g.,  pir,  do-bonding)  molecular  orbitals  at  the  Fermi  energy,  is 
responsible  for  Cooper-Schafroth  pairing  and  superconductivity.  This 

Q 

range  of  3  implies  very  small  values  of  (m/M)  and  p*  [see  Eq.  (60)]  and 

values  of  hwc  approaching  the  Debye  energy  fiwD  [see  Eqs.  (35)  and  (58)]. 

In  the  limit  B  ^  -,  corresponding  to  (m/M)^  (the  "isotope  effect"),  the 
2 

present  model  is  fully  equivalent  ot  the  original  BCS  theory  [see  Eq.  (62)]. 
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5.  For  8  <  1/4,  where  the  relative  displacement  (33)of  the  electron  pair  and 
lattice  ions  is  greater  than  or  equal  to  the  zero-point  vibrational 
amplitude  5q  ^  (m/M)  d  of  the  lattice,  a  static  Jahn-Teller  (or  Peierls) 
distortion  of  the  lattice  can  remove  the  molecular-orbital  degeneracy  at 
the  Fermi  energy,  producing  a  lattice  instability  and  structural 
transformation  (e.g.,  from  cubic  to  tetragonal  symmetry).  An  intrinsically 
low  crystal  symnetry  or  amorphous  structure  can  also  remove  the  molecular- 
orbital  degeneracy  at  the  Fermi  energy.  If  the  energy  gaps  between  the 
resulting  non-degenerate  orbitals  are  not  too  large,  then  Cooper-Scha froth 
pairing  and  superconductivity  can  still  occur  via  dynamic  pseudo-Jahn- 
Teller  coupling  of  the  electrons  and  lattice  ions  induced  by  tne  mixing 

of  the  nearly  degenerate  molecular  orbitals  around  the  Fermi  energy  or 
through  fluctuations  in  orbital  occupancy.  The  above  values  of  3  imply 
relatively  large  values  of  (m/M)13  and  u*  [see  Eq.  (60)]  and  values  of  hwc 
approaching  the  Fermi  energy  [see  Eq.  (58)].  In  the  1 imi t  3  ^  0,  the  present 
model  yields  an  explanation  of  the  vanishing  isotope  effect  in  certain 
(e.g.,  transition-metal)  superconductors,  thereby  precluding  the  necessity 
for  ad  hO£  generalizations  of  BCS  theory  [60-64]  or  esoteric  non-phonon 
mechanisms  of  superconducti vi ty  [10,11,65-67]  in  such  materials.  Finally, 
this  model,  based  on  molecular-orbital  topology  at  the  Fermi  energy  and 
Jahn-Teller  effects,  supports  the  conjecture  that  the  interactions  responsible 
for  lattice  instabilities  and  structural  phase  transformations  are  similar 
to  those  responsible  for  superconductivity  in  the  same  materials  (e.g, 

A15  compounds).  [58] 

6.  the  present  theoretical  model  permits  a  simple  description,  analogous  to 
that  originally  suggested  by  Schafroth,  [56]  of  Cooper  pairs  as 
"quasimolecules"  undergoing  Bose-Einstein-1 ike  condensation  to  the  many- 
electron  superconducting  state.  The  spatially  delocalized  "valence-bond- 
like"  correlation  of  an  electron  pair  over  many  "layered"  or  "tubular" 
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molecular-orbital  components  ip  and  ip_  [see  Figs.  6  and  11  and  Eq.  (6)] 

uses  up  only  a  very  small  fraction  of  the  available  phase  space  defined  by 

the  layered  or  tubular  "cells."  Thus  there  is  ample  phase  space  left 

to  put  more  electron  pairs  into  the  same  valence-bond-like  pair  function  (6). 

As  an  illustration  of  the  quantitative  application  of  the  above  formulae, 

we  continue  with  the  example  of  the  "weakly  Jahn-Teller  coupled"  (8  n#  1) 

2 

superconductor  ,  aluminum,  where  expression  (61)  for  the  superconducting 
transition  temperature  is  valid.  From  the  SCF-Xa-SW  cluster  molecular- 


orbital  contour  map  in  Fig.  6,  which  is  a  direct  measure  of  the  layer  orbital 

O 

parameter  d,  a  value  of  d  ^  1 .8A  appropriate  to  bulk  crystalline  aluminum  is 
obtained.  Substitution  of  the  latter  in  Eq.  (61)  yields  the  transition 


temperature  Tc  ^  1.2°K.  The  latter  agrees  very  well  with  the  experimental 

value  (^c)exp£  x  1-1S°K  for  bulk  aluminum.  [68]  The  Debye  temperature  of  aluminum, 

calculated  on  the  basis  of  expressions  (35)  and  (63),  is  385°K,  which  also 


compares  well  with  the  experimental  value  (Sg)expt  375°K.  [68]  Expression 
(61)  i$  further  consistent  with  the  observed  increase  of  transition  temperature 


with  decreasing  aluminum  particle  size,  namely, 

o  o 


(TcW'5°Kt02'5’Kf0r 


particle  diameters  between  150A  and  50A,  respectively.  [5,6]  Similar  increases 


of  Tc  with  decreasing  grain  size  have  been  measured  for  evaporated  aluminum 
thin  films.  [69,70]  Close  scrutiny  of  Figs.  4  through  6  reveals  a  significant 

O 

increase  of  the  layer-orbi tal  parameter  o  to  a  value  of  ^  2.1A  around  the 


coordinatively  unsaturated  atoms  at  the  cluster  periphery,  due  to  decreased 


net  bonding-anti  bonding  interactions  there,  as  compared  with  the  value 
0 

d  ^  1 .8A  appropriate  to  the  fully  coordinated  central  atom  of  the  cluster 
representing  the  local  environment  of  the  bulk  crystal.  For  particle  diameters 

o  o 

from  150A  to  50A,  the  percentage  of  surface  atoms  with  respect  to  the  total 
number  of  atoms  in  the  particle  varies  from  ^  15%  to  40%,  respectively.  [71] 

Thus,  according  to  expression  (61),  such  increases  of  d  at  small -particle  surfaces 
and  grain  boundaries  should  increase  Tc  to  a  measurable  degree.  Substitution  of 
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the  value  d  ^  2.1A  into  (61)  yields  Tc  ^  1.9°K,  which  is  within  the  range 
(Tc)expt  ^  1-5°K  to  2.5°K  observed  for  particles  between  ^  150A  and  50A.  [5,6] 

O 

Using  weighted  averages  of  the  value  d  'v  2 . 1 A  appropriate  to  atoms  at  the 

0 

cluster  boundary  or  "surface"  and  the  value  d  ^  1 .8A  appropriate  to  bulk 
atoms  in  Eq.  (61)  with  the  weighting  determined  by  the  relative  percentages 
of  surface  and  bulk  atoms,  accounts  for  the  measured  variation  of  T  with 
particle  size.  Note  that,  according  to  Eq.  (35).  such  increases  of  d 
imply  a  reduction  or  "softening"  of  the  phonon  cut-off  frequency  (the 
Debye  frequency)  u>Q,  explaining  the  larger  Tc  within  the  framework  of 
BCS  theory.  [72]  Finally,  comparative  decreases  of  the  orbital  parameter  d 
under  compression  of  the  atoms  seem  to  account  through  Eq.  (61)  for 
substantial  decreases  of  Tc  for  aluminum  under  pressure. 

The  above  theoretical  model  for  superconducting  aluminum  is  also 
consistent  with  the  observation  that  certain  types  of  magnetic  impurities 
in  aluminum  destroy  or  degrade  its  superconductivity.  A  good  example  is 
Mn,  whose  effects  in  aluminum  can  be  modeled  in  a  fashion  analogous  to  the 
treatment  of  Mn  in  copper  (see  Ref.  39)  by  substituting  a  Mn  atom  for  the 
central  A1  atom  in  the  43-atom  cluster  model  of  crystalline  aluminum.  This 
yields  a  MnAl^AlgAl^  cluster  representing  the  local  crystal  environment 
of  the  impurity  up  to  thi rd-nearest  neighbors.  The  SCF-Xa-SW  molecular- 
orbital  energy  levels  for  this  cluster  are  shown  in  Fig.  8,  where  they  are 
compared  with  the  X-ray  photoelectron  spectrum  (XPS)  [73]  of  crystalline 
aluminum  containing  a  small  percentage  of  substitutional  Mn  impurities. 

Comparing  Fig.  8  with  the  SCF-Xa-SW  cluster  moelcular-orbital  energies  for 
pure  aluminum  shown  in  Figs.  2  and  3,  we  see  that  the  main  effects  of  alloying 
Mn. impurities  in  an  aluminum  host,  neglecting  impurity-impurity  interactions,  are 
(1)  to  empty  out  locally  the  otherwise  spatially  extended  tlu(pr)  "layered" 
molecular  orbital  (mapped  in  Figs.  4  through  7)  shown  above  to  be  responsible 
for  the  superconductivity  of  pure  aluminum;  and  (2)  for  the  Mn  3d  orbitals  to 
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overlap  and  hybridize  with  the  surrounding  AT  3p  orbitals  to  produce  a 
narrow  manifold  of  occupied  spatially  localized  bonding  and^  antibonding 
t2g  and  eg  molecular  orbitals  ,  which  are  responsible 

for  the  peak  in  XPS  intensity  just  below  the  Fermi  energy  (the  dashed  line)  in 
Fig.  8.  The  highest  occupied  molecular  orbitals  of  the  MnAl ^ 2A1 gAl 24  cluster 
are  nearly  degenerate  t2^  and  e^  dpo*-dpir*  anti  bonding  orbitals  sharing 
five  electrons,  immediately  above  which  lies  the  emptied  aluminum  tlu(pTr) 
orbital,  as  indicated  in  the  inset  of  Fig.  8.  The  dpo*-dpir*  antibonding  topology 
of  the  former  orbitals  is  clearly  evident  in  the  wave  function  contour  maps 
of  Figs.  9  and  10.  These  orbitals  are  analogous  to  the  t2g  and  eg 

dsdc*-ddsrr*  anti  bonding  orbitals  previously  shown  in  Ref.  39  to  be  associated  with 
the  local  spin  magnetic  moments  of  Mn  and  Fe  impurities  in  a  copper  host  above 
the  Kondo  temperature.  However,  for  Mn  impurities  in  aluminum,  the  fluctuations 
in  occupancy  of  the  nearly  degenerate  individual  t2g  and  eg  antibonding  orbitals 
by  five  electrons,  coupled  with  relatively  small  exchange  splitting  of  the 
orbital  energies,  prevent  stable  magnetic  moment  formation  and  are  responsible 
for  "spin  fluctuations"  in  this  system.  [74]  The  SCF-Xa-SW  cluster  molecular- 
orbital  model  thus  provides  a  real -space  chemical  picture  of  how  Mn  impurities 
in  aluminum  host  degrade  the  superconductivity  of  aluminum  by  locally  emptying 
the  otherwise  spatially  extended  layered  p-r  orbital  responsible  for  Cooper- 
Schafroth  electron  pairing  and  producing  as  the  highest  occupied  orbitals 
spatially  localized  dpo*-dpTr*  antibonding  orbitals  with  flucuating  spin 
magnetic  moment.  These  findings  are  relevant  to  the  general  issue  of  super¬ 
conductivity  vs.  magnetism  which  will  be  discussed  in  Section  VI  of  this  paper, 

IV.  TRAHSITIOU-i-IETAL  SUPERCONDUCTORS 

As  illustrations  of  the  ways  d  orbitals  of  transition  metals,  their 

alloys  and  compounds  can  combine  to  form  around  the  Fermi  energy  coherent, 
spatially  extended  bonding  molecular-orbital  precursors  to  Cooper-Schafroth 


pairing  and  superconductivity,  consider  Figs.  11  and  12.  These  diagrams 
show  schematically  how  d  and  d^2  1  atomic  orbitals,  for  the  chosen  coordinate 
system,  can  overlap  along  a  line  of  atoms  extending  along  the  z-axis  to  form 
"do-type"  molecular  orbitals.  It  is  well  known  in  inorganic  chemistry  [75] 
that  d6  molecular  orbital s  <ontribute  to  multiple  (e.g.  quadruple)  metal-metal  bond 
formation  in  certain  transition-metal  coordination  complexes  (e.g.,  Cl84').  [76] 

In  such  complexes,  the  d<5-bonding  orbitals  are  the  highest  occupied  molecular 
orbitals,  while  the  fully  occupied  dir-  and  dc-bonding  orbitals  have 
significantly  larger  binding  energies.  [77]  Likewise,  d6  orbitals  are  known 
to  be  the  principal  components  of  the  d-band  density  of  states  around  the 
Fermi  energy  in  superconducting  A15  compounds  (e.g.,  Nb^Sn),  [78,79]  where  they 
may  be  viewed  simplistically  as  forming  a  three-dimensional  network  of  inter¬ 
secting  chains  of  d6-bonding  orbitals  like  those  shown  schematically  in 
Figs.  11  and  12.  What  has  not  been  widely  recognized  by  many  solid  state 
physicists  is  that  d5-bonding  molecular  orbitals  of  the  types  shown  schematically 
in  Figs.  13  and  14  are  mainly  responsible  for  the  peak  in  the  d-band  density 
of  states  around  the  Fermi  energy  for  bcc  transition-metal  superconductors 
such  as  V,  Nb,  and  Ta,  as  indicated  in  the  composite  density-of-states  diagram 
of  Fig.  15  based  on  SCF-Xa-SW  cluster  molecular-orbi tal  calculations . 

In  Figs.  11  through  14,  it  can  be  seen  that  the  d5-bond  overlap  is  such 
as  to  form  "tubes"  parallel  to  the  z-axis.  From  these  "tubular  orbitals," 
labeled  <p+  and  according  to  the  sign  of  the  wave  function  d-lobe,  in 
analogy  to  the  p  m  "layer  orbitals"  in  Fig.  6,  one  can  construct  a  Cooper- 
Schafroth  electron-pair  singlet  wave  function  as  a  "valence-bond-like"  product 
of  space  and  spin  components  of  the  form  (6)  and  (7).  Note  in  Figs.  11  through 
14  that,  if  the  tubular  molecular-orbital  components  ;+  correspond  to  a  "spin- 
up"  (t)  electron,  then  the  tubular  components  -i>_  correspond  to  a  "spin-down"  (-) 
electron,  and,  vice  versa,  suggesting  an  "antiferromagnetic"  ordering  of  the 
alternating  di  tubular  orbitals  analogous  to  the  "anti ferromagnetic"  ordering 
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of  prr  layered  orbitals  {fig.  6)  in  superconducting  aluminum  described  in 
Section  III.  However,  the  net  spin  polarization  per  atom  is  zero,  and  there¬ 
fore  this  delocalized  "conduction-electron  anti  ferromagnetism"  should  be 
distinguished  from  conventional  anti  ferromagnetism  localized  on  paramagnetic 
atoms  (see  further  discussion  in  Section  VI).  Indeed,  there  is  recent 
experimental  evidence  for  the  "antiferromagnetic"  ordering  of  conduction 
electrons  in  A15  superconductors.  [80] 

All  of  the  mathematical  steps  in  going  from  the  pair  wave- function  components 
(6)  and  (7)  to  the  final  expression  (56)  for  the  superconducting  transition 
temperature,  derived  on  the  basis  of  p-n  layer-orbital  components  and 
are  also  applicable  to  the  d5  tubular  orbitals  p+  and  <p_.  In  the  latter  case, 
the  interorbital  distance  d  for  the  screened  Coulomb  repulsion  (22)  is  the 
perpendicular  distance  between  parallel  nearest-neighbor  tubular  orbitals 
y+  and  'i>_  at  maximum  d6  bond  overlap.  The  Jahn-Teller  coupling  parameter 
i  is  determined  by  the  do  covalent  overlap  in  much  the  same  fashion  as  the 
p-  covalent  overlap  determines  the  value  of  3  v  j  for  aluminum.  As  already 
mentioned  in  Section  III,  values  of  3  spanning  the  range  0  <  3  <  1/2  are 
appropriate  for  various  transition-metal  systems,  depending  on  d5-bond  over¬ 
lap  and  molecular-orbital  degeneracy  at  the  Fermi  energy. 

While  the  d  (or  degenerate  d  ,  d  )  and  d  2  2  orbitals  in  transition 
xy  3  yz  xz  x  -y 

metal c  can  produce  dc  tubular  molecular  orbitals  of  the  topology  shown  in 
Figs.  11  through  14,  which  are  precursors  to  the  superconducting  state,  2 
orbitals  can  weakly  bond  through  annular  overlap  of  the  type  shown  schematically 
in  Fig.  16.  This  overlap  results  in  "layered"  molecular-orbital  components 
y  and  which,  like  the  layered  p-r  orbitals  shown  in  Fig.  6,  can  be  used 
to.  construct  valence-bond-like  Cooper-Schafroth  electron-pair  wave  functions 
of  the  form  (6).  SCF-Xa-SW  cluster  molecular-orbital  calculations  indicate 
that  overlapping  dz2  orbitals  of  the  type  shown  in  Fig.  16  are  principal 
components  of  the  d  band  at  the  Fermi  energy  in  superconducting  transition- 
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metal  dichalcogenides  (e.g.,  NbSg)  [81]  having  a  layered  crystal  structure. 

As  will  be  shown  below,  somewhat  similar  layered  d^2  orbitals  also  occur  in  the 
vicinity  of  the  Fermi  energy  in  certain  elemental  transition  metals  (e.g.,  Nb) 
and  contribute  along  with  d5-type  orbitals  to  superconductivity. 

A.  Elemental  Transition-Metal  Superconductors  with  bcc  Structure 

With  the  above  general  characteristics  of  do  and  dz2  molecular-orbital 
topologies  in  superconducting  transition-metal  systems  in  mind,  specific 
examples  using  the  formalism  developed  in  Section  III  will  now  be  discussed. 
Consider  first  the  group  of  elemental  transition-metal  superconductors  having 
bcc  crystal  structure,  namely,  V,  Nb,  Ta,  Mo,  and  W.  SCF-Xa-SW  cluster 
molecular-orbital  models  for  these  systems,  along  with  similar  studies  of 
anti  ferromagnetic  bcc  Cr  and  ferromagnetic  Fe,  can  be  summarized  in  the  form 
of  an  approximate  composite  non-spin-polarized  "density-of-states"  profile 
(shown  in  Fig.  15),  in  which  the  locations  of  the  Fermi  energies  of  the 
various  metals  are  indicated  and  the  principal  component  bonding  and 
antibonding  molecular  orbitals  associated  with  peaks  in  the  density  of  states 
are  shown.  Note  in  particular  that  the  Fermi  energies  of  V,  Nb,  and  Ta  are 
located  near  the  density-of-states  peak  associated  with  do-bonding  orbitals  of 
the  types  shown  in  Figs.  13  and  14.  These  Group-VB  elements  are  among  those 
having  the  highest  superconducting  transition  temperatures  in  the  periodic 
table,  namely,  V(TC  *  5.4°K),  Nb(Tc  =  9.2°K),  and  Ta(Tc  =  4.5°K).  [82]  Most 
certainly,  the  high  Tc  values  correlate  roughly  with  the  BCS  formula  (1)  through 
the  nigh  density  of  states  N(0)  at  the  Fermi  energy.  However,  the  variations 
of  Tc  and  significant  deviations  from  the  isotope  effect  among  these  elements 
are  not  consistent  with  the  original  BCS  theory,  [1]  although  formal  extensions 
of. this  theory  have  been  made  to  account  for  these  discrepancies.  [61-64]  As 
will  be  shown  below,  details  of  the  molecular-orbital  topology  at  the  Fermi 
energy,  in  conjuntion  with  expression  (56)  for  T  ,  can  rationalize  the 
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differences  among  the  superconducting  properties  of  V,  Nb,  and  Ta. 

SCF-Xa-SW  cluster  model  calculations  for  vanadium  reveal  that  the  Fermi 
level  coincides  with  a  tri ply-degenerate  ^g^xy’  dyz’  dxz^  molecular  orl3ital 
of  the  do  "tubular"  topology  shown  schematically  in  Fig.  13.  The  actual 
SCF-Xa-SW  numerical  wave  function  contour  map  of  this  orbital  is  displayed 
in  Fig.  17.  From  this  wave  function,  the  intertubular-orbital  distance 

O 

parameter  is  approximately  determined  to  be  d  ^  1.5A,  and  the  Jahn-Teller  coupling 
parameter  appropriate  to  the  d6  bond  overlap  is,  using  expression  (29), 

8  -v  0.25.  Note  that  this  value  of  8  is  on  the  borderline  between  the  ranges 
of  values  appropriate  to  the  dynamic  and  static  Jahn-Teller  effects  [cf.  Eqs.  (26) 
and  (33)].  Indeed,  there  is  some  experimental  evidence  [83]  for  a  cubic-to- 
tetragonal  structural  phase  transformation  (static  Jahn-Teller  distortion)  above 
the  superconducting  transition  to  the  superconducting  state.  Substitution  of  the 
above  values  of  d  and  3  in  expression  (56)  yields  a  value  of  Tc  ^  5.2°K,which 
is  in  good  agreement  with  the  experimental  value  (Tc)exDt  5.4°K. 

For  the  example  of  niobium,  SCF-Xa-SW  cluster  calculations  indicate  that, 
in  contrast  to  vanadium,  the  Fermi  level  coincides  with  a  molecular  orbital 
of  eg(dx2^2,  dz2)  symmetry.  In  V  the  latter  orbital  is  just  below  the 
Fermi  level,  ^g^xy’  dyz’  dxz^’  whereas  in  Nb  this  ordering  is  reversed. 

Contour  maps  of  the  eg(dx2  y2)  orbital  wave  function  at  the  Fermi  energy  of 
Nb  are  shown  for  two  different  planes  in  Figs.  18(a)  and  18(b),  revealing  the 
"tubular"  do  bonding  along  a  chain  of  Nb  atoms.  On  the  other  hand,  the 
degenerate  eg(dz2)  partner  of  this  representation  produces  a  "layered"  orbital 
configuration  of  the  type  shown  in  Fig.  19,  where  the  "pd  hybridization" 
between  nearest-neighbor  Nb  atoms  is  responsible  for  the  complex  topology  of  the 
wave  function.  Clearly,  the  anisotropy  between  the  degenerate  e^(dx2_y2) 

"tubular"  do  orbital  and  eg(dz2)  "planar"  dp-hybrid  orbital  at  the  Fermi 
energy  of  Nb,  as  compared  with  the  topologically  equivalent  degenerate 
t2g(dxy,  dyz>  dxz)  orbitals  of  V,  suggests  that  the  superconducting  properties 


-38- 


of  these  two  transition  metals  are  likely  to  be  intrinsically  different, 
in  spite  of  the  fact  that  they  are  located  in  the  same  column  of  the  periodic 
table  and  have  qualitatively  similar  densities  of  states  (Fig.  15).  Indeed, 
it  has  been  found  by  some  workers  [84]  that  it  is  difficult  to  explain  certain 
tunneling  experiments  on  Nb  without  violating  the  conventions  of  strong¬ 
coupling  superconducting  theory,  whereas  V  seems  to  conform  to  the  latter 

theory.  The  present  theoretical  model  straightforwardly  yields  through  expression 

0 

(56),  for  parameters  d  ^  2.3A  and  Q  ^  0.35  derived  from  the  calculated  SCF-Xa-SW 

molecular  orbitals  of  Nb  and  relation  (29),  a  transition  temperature  of 

Tc  -v  8.6°K,  which  is  in  reasonable  agreement  with  the  experimental  value 

(Tc)expt=  9.2°K.  The  higher  calculated  and  measured  transition  temperatures 

(and  perhaps  other  superconducting  properties)  of  Nb,  as  compared  with  V,  are 

thus  explained  within  the  framework  of  the  present  theoretical  model  as 

resulting  form  a  reversal  in  the  ordering  of  two  nearly  degenerate  topologically 

different  molecular  orbitals,  e^(d^2  ^2,  dz2)  and  t^td^,  d^,  ) ,  at  the 

Fermi  energy,  a  detail  of  the  electronic  structure  which  heretofore  has  not  been 

suspected  on  the  basis  of  conventional  band-structure  density-of-states 

analyses.  [85]  Although  SCF-Xa-SW  molecular-orbital  calculations  have  not  yet 

been  carried  out  for  Ta,  and  therefore  the  exact  ordering  of  e^  and  t^  orbitals 

at  the  Fermi  energy  undetermined  for  this  metal  (relativistic  effects,  e.g., 

spin-orbit  splittings  may  furtner  complicate  this  issue),  it  is  reasonable  to 

expect  on  the  basis  of  crystal  chemistry  that  the  orbital  parameters  3  and 

d  for  Ta  are  more  similar  to  those  for  Nb  than  to  those  for  V.  Assuming  that  the 

Nb  parameters  given  above  are  valid  for  Ta,  one  calculates  from  formula  (56) 

a  value  of  T  ^  2.2°K  if  the  e  (d  2  2,  d  2)  orbital  of  Ta  is  coincident  with 
c  9  x  -y  z 

the  Fermi  energy  and  a  value  of T  ^  6.3°K  if  the  t,„(d  ,  d  d  )  orbital  of  Ta 

L  zg  xy  yz  xz 

is  coincident  with  the  Fermi  energy.  These  calculated  values  bracket  the 
measured  value  (Tc)expt 


=  4.53K  for  Ta. 
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It  will  be  noticed  in  Fig.  15  that  the  Fermi  energies  of  the  Group-VIB 
transition  metals  Cr,  Mo,  and  W  coincide  with  a  region  of  low  density  of  states 
(between  the  bonding  and  antibonding  manifolds)  which  would  normally  be  considered 
to  be  nonbonding.  However,  as  revealed  by  the  contour  map  in  Fig.  20  for  the 
eg(dz2)  cluster  molecular  orbital  wave  function  of  Cr  at  the  Fermi  energy,  this 
low-densi ty-of-states  region  actually  corresponds  to  molecular  orbitals  that  are 
anti  bonding  between  nearest-neighbor  atoms  and  bonding  between  second-nearest 
neighbors.  Since  there  are  eight  nearest  neighbors  and  six  second-nearest  neighbor 
in  the  bcc  structure,  the  net  result  of  this  composite  anti  bonding-bonding  orbital 
topology  at  the  Fermi  energy  is  indeed  mainly  nonbonding.  Similar  molecular- 
orbital  topologies  are  also  found  at  the  Fermi  energies  of  Mo  and  W.  We  suggest 
that  this  type  of  composite  antibonding-bonding  molecular-orbital  topology  at 
the  Fermi  energy  is  responsible  for  the  low-transi tion-temperature  super¬ 
conductivity  of  bcc  Mo  (T  =  0.92°K)  and  W  (T  =  0.01°K),  as  compared  with  bcc 

c  c 

Nb,  through  the  formation  of  spatially  extended  di  "tubular"  orbitals  (cf.  Figs. 

12  and  14),  along  the  direction  perpendicular  to  the  plane  of  Fig.  20,  of 
smaller  "intertubular"  distance  d  than  is  the  case  for  Nb.  These  smaller  values 
of  d  imply  proportional ly  lower  orbital  densities  of  states  at  the  Fermi  energy, 
as  defined  by  Eq.  (49),  and  correspond  to  the  densi ty-of-states  minimum  in 
Fig.  15.  For  Mo,  the  molecular-orbital  parameters  derived  from  the  SCF-Xa-SW 

o  o 

cluster  studies  are  d  a*  1  .3A  and  3  a  0.35  (cf.  d  a  2.3A  and  5  a  0.35  for  Nb), 
which  yield  from  formula  (56)  a  transition  temperature  of  Tc  a  0.3°K  in  good 
agreement  with  the  experimental  value  given  above.  For  W,  preliminary 
relativistic  SCF-Xa-SW  cluster  molecular-orbital  analysis  suggests  that  the 
interorbital  parameter  d  is  further  reduced  bv  the  combined  effects  of  relativity 
and  greater  d-orbital  covalency  along  the  bonding  direction  between  second- 
nearest*neighbor  atoms,  which  accounts  for  the  even  lower  transition  temperature 
of  this  superconductor. 

As  indicated  above,  Cr,  Mo,  and  W  are  all  characterized  by  largely 


-40- 


nonbonding  molecular  orbitals  at  the  Fermi  energy  which  nevertheless  are 
anti  bonding  between  nearest  neighbors  and  bonding  between  second-nearest 
neighbors  (Fig.  20),  yet  Cr  is  an  antiferromagnet  while  Mo  and  W  are  super¬ 
conductors.  It  will  be  discussed  in  Section  VI  of  this  paper  that  for  this  type 
of  composite  molecular-orbital  topology  at  the  Fermi  energy,  the  relative 
amounts  of  spatial  delocalization  of  orbital  charge  between  the  atoms  and 
spatial  localization  of  charge  on  the  atoms  could  largely  determine  whether  the 
metal  is  a  superconductor  or  a  true  antiferromagnet.  Indeed,  it  has  already 
been  argued  above  that  the  spatially  delocalized  correlation  of  Cooper- 
Scha froth  electron  pairs  in  the  tubular  d6  molecular  orbitals  of  the  species 
shown  in  Figs.  11  through  14  corresponds  to  a  type  of  "conduction-electron 
anti  ferromagnetism"  or  "spin-density  wave"  delocalized  between  the  atoms 
(no  net  spin  per  atom),  whereas  conventional  anti  ferromagnetism  is  associated 
with  spin  orbitals  localized  primarily  on  the  atoms  (nonzero  net  spin  per 
atom).  Somewhat  similar  arguments  can  be  used  to  explain  the  occasional 
coexistence  of  superconductivity  and  anti  ferromagnetism  in  the  same  material 
(e.g.,  REMOgSg  Chevrel -phase  compounds  to  be  discussed  in  Section  IVE). 

B.  Elemental  Transition-Metal  Superconductors  with  hep  Structure 

Elemental  transition-metal  superconductors  having  the  hep  structure 

and  an  even  number  of  valence  d  orbitals  per  atom  (e.g.,  Ti ,  Zr,  Ru,  and  Os) 

are  generally  characterized  by  relatively  low  transition  temperatures  (as 

compared,  for  example,  with  bcc  V  and  Nb)  and  vanishing  (or  nearly  vanishing) 

isotope  effects.  [60]  SCF-Xa-SW  cluster  molecular-orbital  studies  of  these 

materials  indicate  that  the  orbital  degeneracy  at  the  Fermi  energy  characterise 

of  cubic  transition  metals  [e.g.,  the  t-  (d  ,  d  ,  d  )  and  e  (d  2,  d  2  2) 

cy  xy  yz  y  z  x  »y 

orbitals  of  V  and  Nb,  respectively,  considered  in  the  preceding  section]  is 
removed  by  the  lower  local  symmetry  of  the  hep  structure,  leaving  fully 
occupied  and  empty  molecular  orbitals,  respectively,  separated  by  a  small  energy 
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gap  around  the  Fermi  energy.  In  other  words,  from  a  real -space  molecular-orbital 
point  of  view,  a  transition  metal  such  as  Zr  can  be  considered  to  have  under¬ 
gone  a  static  Jahn-Teller  distortion  from  the  cubic  to  hep  structure,  removing 
the  orbital  degeneracy  at  the  Fermi  energy,  opening  a  small  energy  gap  there, 
and  consequently  lowering  the  energy  of  the  system.  The  resulting  fully 
occupied  molecular  orbital  just  below  the  Fermi  energy  is  d6-bonding  along 
the  c-axis  between  the  hexagonal  basal  plane  and  triangular  atomic  planes  above 
and  below  the  basal  plane  (as  shown  for  Zr  in  Fig.  21),  whereas  the  resulting 
empty  molecular  orbital  just  above  the  Fermi  energy  is  of  composite  da,  dir- 
bonding  character  within  the  basal  and  triangular  planes  (as  shown  for  Zr  in  Fig. 
22)  and  is  d5-bonding  perpendicular  to  these  planes  (i.e.,  perpendicular  to  the 
plane  of  Fig.  22).  Configurational  mixing  of  these  two  orbitals  will  produce 
"tubular"  do  bonds  around  the  Fermi  energy  like  those  shown  in  Figs.  11  and  12 
spatially  extended  along  the  c-axis  and  passing  continuously  through  the  basal 
and  triangular  planes  of  each  unit  cell.  These  are  the  conditions  described  in 
Section  III  for  Jahn-Teller  coupling  in  the  limit  3  <  1/4  and  for  Cooper- 
Schafroth  pairing  via  dynamic  exchange  of  electrons  between  the  Jahn-Teller-spl it 
molecular  orbitals  immediately  above  and  below  the  Fermi  energy.  For  Zr,  these 
fluctuations  in  orbital  occupancy  around  the  Fermi  energy  correspond  to 

dynamic  displacements  3  of  the  electronic  charge  distribution  at  the  Fermi 

0 

energy  along  the  c-axis  equal  to  5  ^  c/4  -v  1.3A  in  Zr  relative  to  the  average 
lattice  ion  positions.  The  intertubular -orbi tal  distance  parameter  for 
Zr  derived  from  the  calculated  SCF-Xa-SW  cluster  molecular-orbital  wave 

O 

functions  (Fig.  22)  is  d  ^  2.3A.  These  two  quantities,  5  and  d,  are 
consistent  with  Eq.  (33)  only  if  3  ^  0.045,  which  corresponds  to  essentially 
vanishing  isotope  effect.  Substitution  of  these  values  of  2  and  d  into 
formula  (55)  yields  the  transition  temperature  Tc  ^  0.56°K,  which  is  in 
excellent  agreement  with  the  value  (Tc)expt  x  0-55°K  measured  for  pure  Zr. 

It  should  be  noted  that  for  this  very  small  value  of  3,  and  thus  for  the 
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correspondingly  large  value  of  the  quantity  (m/M)  ,  the  repulsive  Coulomb 

interaction  parameter  u*  (or  pseudopotential  U)  defined  in  Eq.  (60)  is 

approximately  60%  of  the  attractive  interaction  parameter  X  (or  pseudopotential 

V)  defined  in  (59).  The  comparative  magnitudes  of  u*  and  X  in  Zr  are 

qualitatively  consistent  with  strong-coupling  extensions  [61-64]  of  BCS  theory 

, « 

which  ascribe  the  vanishing  isotope  effect'  and  low  Tc  for  Zr  (and  similar 
transition  metals)  to  the  effects  of  Coulomb  interactions.  It  may  be  noted 
further  that  for  the  value  of  S  ^  0.045  appropriate  to  Zr,  the  pre-exponential 
"cut-off”  energy  (58)  is  much  larger  than  the  Debye  energy  hu>Q  [see  Eq.  (35)] 
which  corresponds  to  a  value  of  S  ^  1/2.  A  similar  molecular-orbital  analysis 

O  O 

of  hep  Ti  can  be  made,  yielding  the  parameters  5  ^  c/4  'u  1.2A,  d  'v  2.1A, 
and  6  ^  0.050,  again  consistent  with  the  vanishing  isotope  effect.  Substitution 
of  the  last  two  quantities  into  (56)  gives  Tc  -v  0.39°K,  in  (clearly  fortu’ tously) 
exact  agreement  with  the  measured  value  (Tc)expt  =  0.39°K. 

The  same  theoretical  approach  can  also  be  applied  to  the  hep  transition- 
metal  superconductors  Ru  and  0s,  although  because  of  their  positions  further 
toward  the  right  of  the  periodic  table  and  hence  their  greater  number  of 
valence  d  orbitals  per  atom,  the  do-bonding  molecular  orbitals  around  the 
Fermi  energy  responsible  for  superconductivity  have  a  somewhat  different 
orientation  from  those  described  above  for  Zr  and  Ti.  Nevertheless,  the 

O 

molecular-orbital  parameters,  d  ^  2.0A  and  3  'v  0.05,  calculated  for  both 
Ru  and  0s  are  of  the  same  order  of  magnitude  as  those  given  above  for  Zr 
and  Ti.  Substitution  of  these  parameters  in  Eq.  (56)  yields  the  transition 
temperatures,  Ru(T  'u  0.43°K)  and  0s(T  ^  0.64"K),  consistent  with  the 
respective  measured  values  Ru[(Tc)expt  *  0.49°K]  and  0s[(Tc)expt  =  0.66°K]. 

The  only  major  descrepancy  between  theory  and  experiment  here  is  that, 
while  both  indicate  effectively  vanishing  isotope  effect  for  Ru,  the  measured 
dependence  of  Tc  on  atomic  mass  corresponds  to  a  value  of  (2)  t  s  0.20, 
whereas  the  calculated  value  is  U)theor  ^  0.05.  However,  since  ail  four 
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hcp  transition  metals,  Ti,  Zr,  Ru,  and  Os,  have  similar  bonding  configurations, 

transition  temperatures,  and  Jahn-Teller  coupling  parameters  B,  it  is  possible 

that  the  observed  larger  isotope  effect  in  Os  is  an  artifact  of  materials 

preparation.  In  any  case,  the  all-over  agreement  between  the  present 

theoretical  model  and  experiment  for  these  superconductors  is  somewhat 

better  than  that  obtained  from  "strong-coupling"  BCS-type  theories.  [61-64] 

The  hep  transition-metal  superconductors,  Tc  and  Re,  deserve  separate 

discussion  from  Ti,  Zr,  Ru,  and  Os,  because  they  have  comparatively  high 

transition  temperatures,  Tc(T  =  7.3°K)  and  Re(T  =  1.7°K),  and  in  the  case 

c  c 

of  Re  where  isotope  measurements  have  been  carried  out,  a  comparatively 
large  isotope  effect  corresponding  to  (S)expt  =  0*39.  [64]  These  differences 
are  explained,  first  of  all,  by  the  fact  that  Tc  and  Re  have  an  odd,  rather 
than  even,  number  of  valence  d  electrons  per  atom,  so  that  the  d6  molecular- 
orbital  degeneracy  at  the  Fermi  energy  is  not  completely  removed,  despite  the 
non-cubic  structure.  Consequently,  these  materials  are  subject  to  ordinary 
dynamic  Jahn-Teller  coupling  in  the  range  of  1/4  <  3  <  l/2[see  Eq.  (26)]. 

v  J 

On  the  basis  of  SCF-Xa-SW  cluster  molecular-orbital  analysis,  one  obtains  for 

0 

Tc  values  of  3  'u  0.30  and  d  ^  1.9A,  yielding  from  formula  (56)  a  transition 
temperature  of  Tc  8.1°K.  For  Re,  the  calculated  molecular-orbital  parameters 

O 

are  3  %  0.40  and  d  n,  1.9A,  giving  Tc  ^  2.1°K.  Both  theoretical  values  of  the 
transition  temperature  are  in  good  agreement  with  the  measured  values  for 
Tc  and  Re  given  above,  and  the  calculated  values  of  B  are  consistent  with 
the  measured  isotope  effect. 

In  conclusion  to  this  and  the  preceding  section  on  elemental  transition- 
metal  superconductors,  it  should  be  pointed  out  that  the  well-known  correlation 
of  Tc  with  electron-to-atom  ratio  ,  originally  emphasized  by  Matthias  [3]  and 
discussed  by  many  other  workers  [12-15,60-64]  can  be  understood  in  terms 
of  the  molecular-orbital  character,  occupancy,  and  degeneracy  at  the  Fermi 
energy,  as  described  above. 
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C.  A15  Transition-Metal  Compound  Superconductors 

Transition-metal  compounds  such  as  Nb^Sn  and  Nb^Ge  crystallize  in  the 
A15  structure,  where  the  Nb  atoms  form  a  three-dimensional  network  of 
orthogonal  linear  chains  extending  along  the  faces  of  the  cubic  unit  cells 
and  the  Sn  or  Ge  atoms  are  located  at  the  corners  and  center  of  each  unit 
cell.  These  compounds  are  among  the  highest  temperature  superconductors 
known,  e.g.  NbjSntT.  *  18.0°K)  and  Nb^Ge^  «  23.26K).[86]  SCF-Xa-SW  cluster 
molecular-orbital  studies  of  Nb^Sn  reveal  the  presence  of  spatially 
delocalized  "tubular"  dd-bonding  molecular  orbitals  In  the  vicinity  of 
the  Fermi  energy  of  the  types  shown  schematically  in  Figures  11  and  12 
and  quantitatively  in  Figure  23  along  four  atoms  of  each  chain.  While 
these  tubular  d6  orbitals  at  the  Fermi  energy  will  be  shown  below  to  be 
solely  responsible  for  the  superconductivity  of  Nb^Sn  (and  other  A15 
compounds),  they  are  not  the  sole  contributors  metal -metal  bonding  in 
these  materials,  which  also  depends  on  do  and  dr-bonding  molecular  orbitals 
that  are  spatially  localized  between  pairs  of  Nb  atoms  in  each  chain  and  lie  well 

below  the  Fermi  energy  as  they  do  in  elemental  bcc  crystalline  Nb  (see  Fig.  15). 

Indeed,  the  effective  electron  configuration  of  a  pair  of  nearest-neighbor  Nb 

2  4  2 

atoms  along  each  chain  in  Nb^Sn  is  close  to  that,  (dc)  (dr)  (d5)  ,  for 

4- 

"quadruply-bonded"  coordination  complexes  (e.g.  Mo2Clg  ^  ’nor9anic 
chemistry.  [75-77]  Although  the  Sn  atoms  and  neighboring  Nb  chains  are  not 
explicitly  shown  in  Figures  11,  12,  and  23,  the  d<5  bonding  (as  well  as  the 
do  and  dr  bonding)  along  each  chain  is  due  not  only  to  the  direct  covalent 
overlap  between  nearest-neighbor  Nb  4d  orbitals  but  is  also  promoted  by  the 
interaction  and  hybridization  between  Nb  4d  orbitals  and  neighboring  Sn  5p 
orbitals  and  between  the  4d  orbitals  of  Nb  atoms  on  neighboring  chains.  The 
net  effect  of  these  interactions  is  to  produce  a  Nb-Nb  nearest-neighbor  bond 

O 

distance  (2.64A)  along  each  chain  in  NbgSn  which  is  significantly  less  than  the 
nearest-neighbor  distance  (2.86A)  characteristic  of  elemental  bcc  crystalline 
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Nb.  The  influence  of  the  Sn  atoms  on  Nb-Nb  bonding  is  analogous  to  that  of  the 

4- 

C1  ligands  on  Mo-Mo  bonding  in  the  molecular  coordination  complex  MOgClg  , 

0 

where  the  Mo-Mo  bond  distance  (2.14A)  is  even  shorter  than  the  nearest- 

0 

neighbor  distance  (2.73A)  characteristic  of  elemental  bcc  crystalline  Mo. 

As  will  be  shown  below,  this  type  of  ligand-metal  interaction  is  a  key 
contributing  factor  in  understanding  the  variation  of  the  transition  temperature, 
which  is  related  to  the  d6  bonding  at  the  Fermi  energy  through  parameters 
d  and  3  in  (56),  among  the  A15  superconductors. [86]  Thus  it  should  be 
emphasized  that  Figs.  11,  12,  and  23  do  not  imply  uncoupled  one-dimensional 
d6-bonded  Nb  chains  as  in  the  simple  qualitative  Labbe-Friedel  [87]  band- 
structure  model  for  A15  compounds,  but  instead  a  three-dimensional  network 
of  coupled  df-bonded  Nb  chains  including  the  effects  of  Nb(4d)-Sn(5p)  metal- 
ligand  interactions,  consistent  with  the  results  of  more  recent  quantitative 
band-structure  calculations.  [78,79,88]  The  predicted  spatial  delocalization 
of  metal -metal  bond  charge  along  the  Nb  chains  in  Nb^Sn  is  qualitatively 
consistent  with  x-ray  data  [89]  showing  an  excess  of  charge  along  the  metal  chains 
in  the  A15  compound  V^Si,  although  in  the  latter  compound  somewhat  less  d<5-bond 
overlap  is  predicted  by  the  present  theoretical  model  because  of  the  greater 
spatial  localization  of  V  3d  orbitals  in  comparison  to  Nb  4d  orbitals. 

As  in  the  preceding  examples  of  elemental  transi tion-metal  superconductors, 
the  three-dimensional  networks  of  tubular  do  molecular  orbitals  at  the  Fermi 
energy  of  the  types  shown  in  Figs.  11  and  12  can  be  used  to  construct  a 
Cooper-Schafroth  singlet  pair  function  as  a  "valence-bond-like"  product  of 
space  and  spin  components  of  the  form  (6)  and  (7).  Note  in  Figs.  11  and  12  the 
spatially  delocalized  "anti ferromagnetic"  ordering  of  electron-pair  spins  along  the 
"tubular"  do  bonds,  with  a  net  spin  per  atom  equal  to  zero,  as  discussed  in 
preceding  sections.  Indeed,  there  is  recent  experimental  evidence,  in  the  form 
of  conduction-electron-spin-resonance  (CESR)  measurements,  for  "anti ferromagnetic" 
ordering  of  conduction  electrons  in  the  A15  compound,  Nb^Ge.  [80]  This  type 
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of  "anti ferromagnetic"  ordering  of  conduction  electrons  at  the  Fermi  energy 
should  not  be  confused  with  the  conventional  anti  ferromagnetism  of  spatially 
localized,  paramagnetic  d  electrons,  nor  should  it  be  confused  with  the 
"coexistence"  of  anti  ferromagnetism  and  superconducti vity  observed  in  certain 
other  transition-metal  compounds  (see  Section  I V  E) .  The  anti  ferromagnetic 
ordering  in  A15  superconductors  can,  in  fact,  be  viewed  as  a  conduction- 
electron  spin-density  wave  of  the  type  predicted  by  Overhauser.  [90] 

It  has  been  customary  to  calculate  Tc  values  for  A15  compounds  using  the 
BCS  formula  (1)  and  the  total  electronic  density  of  states  at  the  Fermi 
energy,  N(0)  derived  from  band-structure  calculations.  [91,78,79]  The 
occupied  d-band  densities  of  states  of  these  compounds  can  be  partitioned 
into  da-,  dir-,  and  d6-bonding  molecular-orbital  components,  as  for  the 
example  of  bcc  transition  metals  in  Fig.  15,  with  the  do*-,  dir*,  and  dc*- 
antibonding  components  lying  above  the  Fermi  energy.  [78]  The  Fermi  energies 
themselves  generally  coincide  with  a  density-of-states  peak  corresponding 
primarily  to  a  three-dimensional  array  of  intersection  chains  of  "tubular" 
dc(d  2  2,  d  2  2,  d  2  2)  orbitals  or  the  type  shown  in  Fig.  12.  For  Nb^Sn, 

A*  Jr  Y  "  fc  A  "  4  v 

the  do-bond  overlap  at  the  Fermi  energy,  as  revealed  by  SCF-Xa-SW  studies, 
is  substantial,  implying  frcm  expression  (29)  a  value  of  5  ^  0.15. 

O 

For  the  derived  interorbital  distance  d  x  1.9A  [the  maximum  distance  between 
neighboring  tubular  functions  and  in  Fig.  12]  in  Nb-^Sn,  formula  (56)  yields 
Tc  -v  18.4°K.  The  relatively  small  value  of  the  Jahn-Teller  coupling  parameter 
3  is  consistent  with  the  nearby  vanishing  isotope  effect  in  Nb^Sn  and  further 
suggests  that  the  static  Jahn-Teller  effect  can  compete  with  Cooper-Schafroth 
electron  pairing,  resulting  in  a  structural  phase  transformation  of  the  type 
mentioned  in  Section  III.  Indeed,  a  distortion  of  the  lattice  from  cubic  to 
tetragonal  structure  is  known  to  occur  in  Nb^Sn  and  other  high-Tc  AI5  super¬ 
conductors  and  is  suspected  to  lower  the  superconducting  transition  temperature 
from  that  otherwise  attainable  for  the  cubic  structure.  [58]  A  reduction  of 
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the  intertubular  orbital  distance  d  in  Eq.  (56)  with  such  a  lattice  distortion 
would  be  consistent  with  a  reduction  in  the  value  of  T  . 

D.  "Layered-Structure11  Trans tion-Metal  Compounds 

Another  interesting  class  of  superconducting  transition-metal  compounds 
are  the  layered  transition-metal  dichalcogfindides,  such  as  Nb$2»  Ta$2»  NbSe2. 
and  TaSej.  [81]  SCF-Xa-SW  cluster  molecular-orbital  calculations  for  these 
compounds  indicate  that  the  Fermi  energy  is  coincident  with  molecular  orbitals 
of  d^2  symmetry,  shown  schematically  in  Fig.  16  for  three  neighboring  metal 
atoms  in  a  single  layer.  This  orbital  is  moderately  anti  bonding  between  the 
transition-metal  d  orbitals  and  chalcogen  "ligand"  p  orbitals  (not  shown  in 
Fig.  16)  but  weakly  bonding  between  the  "doughnut-like  annular  regions  of  the 
metal  dz2  orbitals  in  each  layer  (the  dashed  profiles  shown  in  Fig.  16).  The  bond¬ 
ing  between  annular  regions  forms  a  diffuse  but  spatially  coherent  wave  function 
throughout  each  layer,  thereby  satisfying  the  orbital  criterion  for  the  existence 
of  superconductivity.  In  this  respect,  the  superconducting  state  of  the  layered 
transition-metal  dichalcogenides  is  "two-dimensional."  However,  the  total 
superconducting  state  of  the  crystal  is  a  composite  result  of  the  three- 
dimensional  array  of  layered  dz2  orbitals,  just  as  the  superconducting  states 
of  aluminum  and  A15  compounds  discussed  above  are  the  composite  results  of 
three-dimensional  networks  of  layered  ot  and  tubular  d5  orbital  con¬ 
figurations.  The  relatively  low  transition  temperatures  of  these  layered 
superconductors,  e.g.,  NbS2(T.  *  68K)  and  TaSe2  (T£  *  0.2°K),  as  compared  with 
those  of  the  d<$-bonded  AT 5  compounds,  e.g.,  Nb-jSn  (T£  *  18°K)  and  V^Sn  (Tc  *  3.8°K) 
are  consistent  with  the  fact  that  the  bond  overlap  between  the  annular  parts 
of  the  dz2  orbitals  in  the  former  compounds  (see  Fig.  16)  is  significantly  less 
than  that  between  the  d6  orbitals  in  the  latter  compounds  (cf.  Fig.  11).  The 
exact  amounts  of  bond  overlap  and  concomitant  Tc  values  for  the  layered  super¬ 
conductors  are  dependent  on  the  anti  bonding  contributions  of  the  chalcogen 


ligand  p  orbitals  to  the  dz2  molecular  orbitals,  in  much  the  same  fashion  as 
the  p-r  bonding  orbital  components  in  aluminum  are  promoted  by  the  pc*  anti- 
bonding  components  (see  figs.  4  through  6),  and  in  much  the  same  way  as  the  net 
d6  bonding  along  the  chains  of  metal  atoms  in  the  A15  compounds,  Nb3$n  and 
V^Sn,  is  dependent  on  the  interaction  of  the  metal  atoms  with  the  neighboring 
Sn  ligands.  The  intercalation  of  organic  molecules,  such  as  pyridine,  between 
the  layers  of  superconducting  transition-metal  dichalcogenides  does  not 
appreciably  alter  their  Tc  values.  [81]  This  is  explained  by  the  fact  that 
the  bonding  between  the  layers  and  intercalate  molecules  occurs  via  the  dz2 
lobes  perpendicular  to  each  layer  (the  solid  orbital  contours  in  Fig.  16),  which 
leaves  largely  unaffected  the  bonding  between  annular  regions  within  each 
plane  (the  dashed  orbital  contours  in  Fig.  16)  responsible  for  the  super¬ 
conducting  state.  The  weak  annular  type  of  bonding  between  dz2  orbitals  is 
not  strictly  limited  to  layered  transition-metal  compounds  but  can  also 
contribute,  along  with  the  more  dominant  dc  bonding,  to  the  superconducting 
state  in  transition  metals  and  their  alloys .(see  Section  IVA). 

E.  "Chevrel -Phase"  Molybdenum  Sulfide  Compounds 

Superconducting  molybdenum  sulfides  having  the  composition  RExMo6Sg 
(Chevrel  compounds)  containing  a  high  concentration  (x  ^  1)  of  magnetic 
(RE)  ions  [92]  are  of  considerable  interest  because  of  the  apparent  coexistence 
of  superconductivity  and  magnetism.  [93]  A  preliminary  molecular-orbital 
analysis  of  the  MOgSg  octahedral  cluster  which  is  the  key  structural  unit  in 
the  superconducting  Chevrel  phases  has  been  carried  out,  using  the  results 
of  recent  SCF-Xa-SW  studies  of  similar  Mo-atom  cluster  compounds  by  Cotton  and 
Stanley.  [94]  This  analysis  suggests  that  the  superconducting  state  of  these 
Chdvrel  phases  and  the  more  common  ones  not  containing  Re  ions  [95]  arises 
exclusively  from  coherent  spatially  delocalized  d-orbital  bonding  between 
the  Mo  atoms,  promoted  in  part  by  the  direct  covalent  overlap  between  Mo  d 
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orbitals  in  the  Mo^Sg  clusters  and  in  part  by  the  antibonding  contribution  of 
the  interaction  between  Mo  d  orbitals  and  S  ligand  p  orbitals.  The  latter 
contribution  may  be  compared  with  the  above-described  effect  of  the  pc*- 
antibonding  component  of  aluminum  in  promoting  coherent  spatially  extended 
ptr  bonding  (see  Figs.  4  through  6).  On  the  other  hand,  the  magnetic  properties 
of  the  Re^MOgSg  compounds  are  associated  exclusively  with  localized  f  orbitals 
on  the  RE  ions  which  are  formally  but  weakly  anti  bonding  with  respect  ot  the 
MOgSg  clusters.  Thus  the  interaction  between  the  magnetic  RE  ions  and  the 
superconducting  electrons  is  weak,  explaining  why  the  presence  of  such  ions  does 
not  generally  destroy  the  superconducti vity  of  the  Chevrel  phases. 

F.  Palladium  Hydride  (PdH^) 

Pure  palladium  is  strongly  exchange  enhanced  but  not  superconducting.  [64] 

On  the  other  hand,  the  addition  of  large  amounts  of  interstitial  hydrogen  to 
palladium  to  form  the  alloy,  palladium  hydride  (PdHx),  eliminates  the  spin 
fluctuations  and  produces  a  good  superconductor  (e.g.  2.1 3K  <  T  <  9.1°K  for 

'  u  C  ^ 

0.35  <  x  <  1 ) .  [96]  SCF-Xa-SW  molecular-orbi tal  studies  of  a  19-atom 

palladium  cluster  representing  the  local  molecular  environment  of  pure 

crystalline  palladium  up  to  second-nearest  neighbors  indicate  the  presence 

of  degenerate  spin-polarized  dd~*  and  dpc*  antibonding  orbitals  at  the  Fermi 

energy.  [97]  Fluctuations  of  electron  occupancy  between  these  two  types  of 

degenerate  orbitals  are  discrete  cluster  analogs  of  the  observed  magnetic 

spin  fluctuations.  Similar  SCF-Xa-SW  cluster  molecular-orbi tal  studies  of  PdH  [9 

x 

indicate  that  the  interstitial  hydrogen  in  palladium  raises  the  Fermi  energy 
above  the  Pd-Pd  ddr*  and  dpo*  anti  bonding  orbitals  and  makes  it  coincident  with 
orbitals  having  spatially  delocalized  components  which  are  so-bonding  within 
and  so*-antibonding  between  parallel  planes  of  weakly  interacting  H  atoms 
(see  Fig.  24),  thereby  quenching  the  spin  fluctuations.  The  present  theoretical 
model  further  suggests  that  these  so-bonded  planes  of  H  atoms  form  layer 
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orbitals  <|/+  and  i at  the  Fermi  energy  (Fig.  24)  which  are  somewhat  similar 
to  those  for  aluminum  (cf.  Fig.  6),  and  that  the  valence-bond-like  correlation  of 
electrons  in  these  layer  orbitals  [see  Eq.  (6)],  in  conjunction  with  with  the 
dynamic  Jahn-Teller  effect  is  responsible  for  Cooper-Schafroth  electron  pairing 
and  superconductivity  in  PdHx- 

V .  OTHER  SUPERCONDUCTORS 

A.  Polymeric  and  Organic  Superconductors 

Prime  examples  of  superconducting  polymeric  and  organic  solids  are 
polysulfur  nitride  ( SN ) x  and  tetramethyltetraselenafulvalene-hexafluorophosphate 
(TMTSF^PF^..  Crystalline  (SN)x  is  not  only  a  good  ordinary  electrical  conductor 
along  the  direction  of  "fiber  bundles"  formed  by  large  aggregates  of  individual 
polymer  molecules,  but  is  also  an  anisotropic  superconductor,  albeit  a  low- 
T  0.25°K)  one. [99]  Previous  SCF-Xct-SW  studies  [100]  reveal  that  the 
Fermi  energy  is  coincident  with  coherent  p--bonding  molecular  orbitals  spatially 
delocalized  along  the  ( SN ) x  chains.  Despite  the  one-dimensional  nature  of  the 
individual  chains,  they  result  in  an  effective  three-dimensional  array  of 
p-r-bonding  orbitals  due  to  the  pT  interchain  coupling  in  the  (100)  plane  and 
pc  interchain  coupling  in  the  (T02)  plane.  [100]  These  orbitals  are  somewhat 
analogous  to  the  composite  p~-pc*  molecular  orbitals  (Figs.  4  and  5)  shown  in 
Section  III  to  be  precusors  to  superconductivity  in  aluminum.  However,  in  ( SN ) x 
the  p^-bonding  components  are  exclusively  along  the  direction  of  the  fibers  and 
thus  lead  to  anisotropic  superconductivity,  wheraas  in  aluminum  there  are  three 
principal  cubic  crystallographic  directions  which  result  in  isotropic  superconduc¬ 
tivity.  The  anisotropic  prr  bonding  and  relatively  weak  pr  bond  overlap  associated 
with  the  "zig-zag"  structure  of  the  (SN)x  polymers  are  likely  responsible  for  the 
low  T£  value  of  this  material. 

The  organic  charge-transfer  salt,  (TflTSF^PFg,  the  crystal  structure  of 
which  consists  of  stacks  of  TMTSF  molecules  ordered  in  sheets  separated  by  PFg 
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sheets  (see  Fig.  25),  is  observed  to  be  superconducting  at  Tc  ^  0.9°K  for 
an  applied  presure  of  12k  bar.  [101]  A  structurally  similar  compound, 

(TMTSF) gCl 04 ,  has  more  recently  been  found  to  be  superconducting  at  Tc  ■v  1.2°K 
without  any  applied  pressure[102]  A  preliminary  molecular-orbital  analysis, 
based  on  previous  SCF-Xa-SW  calculations  [103]  for  TTF  dimers  and  trimers 
in  the  somewnat  similar  (but  nonsuperconducting)  compound,  TTF-TCNQ,  suggests 
that  Coo  per -Sc  ha  froth  pairing  and  superconductivity  in  (TMTSF) 2PFg  and 
(TWTSF)2C104  arise  from  the  "valence-bond-1 ike"  pairing  [see  Eq.  (6)]  of 
pi-bonding  "layer-orbital"  components  ip  and  at  the  Fermi  energy  analogous 
to  those  shown  in  Figs.  4  through  7  for  aluminum.  In  (TMTSF)2PFg  and 
(TMTSF)2C104,  the  layer  orbitals  and  originate  from  the  combined  effects 
of  interacting  pi-bonded  TMTSF  molecules  in  each  stack  and  "phase-locked" 
coupling  between  TMTSF  stacks  via  weak  overlap  with  the  otherwise  "nonbonding" 
pi  orbitals  localized  on  the  F  and  0  ligands,  respectively,  of  the  intervening 
PFg"  and  Cl 0 4 _  anions.  The  layer-orbital  components  and  are  shown 
schematical ly  for  (TMTSF)2PFg  in  Fig.  25,  including  the  "anti ferromagnetic" 
correlation  of  conduction-electron  spins  between  the  molecules,  in  analogy 
to  the  spin  correlations  of  the  layer  orbitals  of  aluminum  shown  in  Fig.  6. 
Indeed,  there  is  some  experimental  evidence  for  such  anti  ferromagnetic  ordering 
in  (TMTSF)2PFg,  corresponding  to  spin-density  waves.  [104]  It  is  further 
suggested  that  in  the  crystal,  (TMTSF)2PFg,  the  TMTSF-PFg-TMTSF  pi-bond 
overlap  necessary  for  phase-locked  interchain  coupling  and  hence  for  Cooper- 
Schafroth  pairing  of  spatially  extended  layer  orbitals  <{/+  and  cannot 
be  attained  unless  external  pressure  reducing  the  TMTSF-PFg  intermolecular 
distance  is  applied,  whereas  in  (TMTSF)2C10^  these  conditions  are  present 
without  applied  pressure.  Under  these  conditions,  the  interlayer  (<*/+-  <p_) 

O 

orbital  distance  in  these  components  is  d  ^  1.8A,  i.e.  approximately  half 
the  TMTSF  intermolecular  distance  in  a  stack  (see  Fig.  25).  The  pi-bond 
overlap  is  small,  so  that  according  to  expression  (29)  the  Jahn-Teller  coupling 
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parameter  i s  3  ^  j.i-e-the  limit  of  weak  dynamic  Jahn-Teller  (electron-phonon) 
coupling.  Substitution  of  these  parameters  into  formula  (56),  using  for  M  weighted 
averages  of  the  masses  of  the  atoms,  [Se,  C,  and  F  in  (TMTSF)2PFg]  and  [Se,  C,  and 
1  in  (TMTSF^ClOjj],  participating  in  layered  pm  bonding  yields  0.8°K  <  Tc  <  1.2°K, 
which  is  in  excel >ent  agreement  with  the  range  of  measured  values  (Tc)expt  ^  0.9°K 
and  1.2°K  for  (TMTSF)2PFg  and  (TMTSF)2C104,  respectively.  The  absence  of 
superconductivity  and  presence  of  charge-density  (instead  of  spin-density)  waves 
in  TTF-TCNQ  [104]  is  most  likely  due  to  the  "herringbone"  structure  (Fig. 26)  of  the 
TTF  and  TCNQ  chains,  which  prevents  the  two-dimensional  interchain  coupling 
and  pahse  locking  of  the  pm  electrons  and  thus  does  not  meet  the  conditions 
[present  in  (TMTSF)2PFg]  for  spatially  extended  layer  orbital  formation. 

B.  Lanthanides  and  Actinides 

Of  the  lanthanides,  only  La  is  superconducting,  while  there  are  three 
actinides,  Th,  Pa,  and  V  which  are  superconducting.  [64]  Of  these,  the 
superconductivity  of  u-uranium  is  most  anomalous  because  of  the  sensitivity  of 
Tc  to  pressure  and  because  of  the  large,  positive  isotope  effect.  [105]  The 
latter  effect  has  been  presented  as  evidence  for  a  "nonphonon"  mechanism  of 
superconductivity  in  x-U  involving  virtual  excetations  of  electrons  to  low- 
lying  f  states.  [105]  On  the  other  hand,  theoretical  arguments  have  been 
presented  to  show  that  the  large,  positive  isotope  effect  in  a-u  can  be 
explained  within  the  framework  of  BCS  theory.  [106]  Preliminary  SCF-Xa-SW 
cluster  molecular-orbital  models  for  a-ll  indicate  that  the  superconductivity 
of  this  material  is  mainly  associated  with  do-bonding  orbitals  and  dynamic 
Jahn-Teller  instability  at  the  Fermi  energy,  but  with  important  effects  due 
to  'mixing"  with  or  excitations  (of  the  type  described  in  Section  III)  to 
nearly  degenerate  f  orbitals.  Thus  these  preliminary  results  support  the 
original  conjecture  [105]  of  virtual  excitations  to  low-lying  f  orbitals 
in  a-u,  although  it  is  the  d<5-orbital  topology  at  the  Fermi  energy  and  associated 
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dynamic  Oahn-Teller  electron-lattice  ("phonon")  coupling  that  is  the 
mechanism  for  superconductivity  in  a-U,  consistent  with  the  conclusions  of 
Ref.  106. 

C.  Amorphous  Alloys 

Since  the  first  superconducting  amorphous  alloy  obtained  by  rapid  liquid 
quenching  was  reported  in  1975,  a  systematic  search  for  new  superconducting 
amorphous  alloy  compositions  has  been  made.  Approximately  fifteen  amorphous 
alloys  are  now  known  to  be  superconductors  with  transition  temperatures  ranging 
from  1.5  to  9.0°K.  [107]  Among  these  are  zirconium -based  metallic  glasses 
such  as  Zr65Cu35(Tc  =  2.0°K),  Z^Ni^T.  =  2.9“K),  and  Zr70Co3()(Tc  *  3.3°K).  [108] 
Of  course,  there  are  many  amorphous  alloys,  such  as  the  Si-Pd  system,  which 
are  not  superconducting. 

Recently,  SCF-Xa-SW  cluster  molecular-orbital  models  for  the  local 
electronic  structures  of  Zr-Cu  and  Si-Pd  metallic  glasses  have  been 
constructed.  [109]  Both  alloys  are  characterized  by  a  fully  occupied  d  band 
associated  primarily  with  the  Cu  or  Pd  components  and  a  partially  occupied 
Zr  d  band  or  Si  p  band  lying  at  higher  energies,  although  there  is  significant 
Cu-Zr  or  Pd-Si  bonding  and  Zr-Cu  or  Si-Pd  antibonding  character  within  these 
respective  bands  and  some  electron  transfer  from  Zr  to  Cu  or  from  Si  to  Pd. 

The  calculated  electronic  structures  are  in  good  quantitative  agreement  with  and 
provide  an  interpretation  of  published  photoelectron  spectra  for  these  alloys. 

Perhaps  the  most  striking  feature  of  the  calculated  electronic  structures  is 
the  nature  of  the  electronic  states  coincident  with  or  just  below  the  Fermi 
energy,  i.e.  the  highest  occupied  molecular  orbitals.  These  states  can  be 
described  as  forming  highly  delocalized  and  spatially  directed  Zr-Zr  d-orbital 
or  Si -Si  p-orbltal  chemical  bonds  which  are  promoted,  in  part,  by  the 
composite  bonding-anti  bonding  interaction  with  the  d  orbitals  localized  on 
neighboring  Cu  or  Pd  components.  These  bonds  are  clearly  visible  in  the 


molecular-orbital  wave-function  contour  maps  shown  in  Fig.  27  for  the  Zr-Cu 
and  Si -Pd  systems,  respectively. 

The  existence  of  directed  Zr-Zr  and  Si -Si  bonds  in  amorphous  Zr-Cu  and 
Si-Pd  alloys,  respectively,  may,  at  first  right,  seem  surprising  since  the 
Zr  and  Si  atoms  are  not  nearest  neighbors  itf  these  alloys  and  the  Zr-Zr  and  Si- 
Si  bond  distances  are  significantly  larger  than  those  in  elemental  crystalline 
Zr  and  Si.  Nevertheless,  the  effects  of  the  nearest-neighbor  Cu  and  Pd 
d  orbitals  in  acting  as  "bridges"  for  the  formation  of  Zr-Zr  and  Si -Si  bonds  are 
clearly  evident  in  Fig.  27. 

On  a  scale  spanning  many  atoms,  the  Zr-Zr  and  Si -Si  bonds  at  the  Fermi 
energy  form  a  topologically  disordered  "network"  which  is  believed  to  be 
related  to  the  glass-forming  capabilities  of  these  and  similar  alloys.  [110] 

In  the  case  of  amorphous  Zr-Cu  alloys,  the  Zr-Zr  bonds  at  the  Fermi  energy, 
despite  the  lack  of  long-range  crystalline  order,  are  of  the  do  type  shown 
in  Section  IV  to  be  precursors  to  superconductivity  in  transition  metals 
and  their  compounds.  On  the  other  hand,  the  Si-Si  bonds  at  the  Fermi  energy 
of  amorphous  Si-Pd  alloys  are  of  the  more  localized  pa  type  which  is  not 
favorable  for  Cooper-Schafroth  pairing  and  superconductivity.  Indeed,  amorphous 
Zr-Cu  alloys  are  superconductors,  [108]  whereas  amorphous  Si-Pd  alloys  are  not. 

VI.  MAGNETISM  VS.  SUPERCONDUCTIVITY 

Although  magnetism  and  superconductivity  can  occasionally  coexist  in 
some  complex  materials  (see  Section  IV  E),  it  is  more  generally  true  that 
superconductivity  and  magnetism  are  mutually  exclusive  in  a  material  at  the 
same  temperature  and  pressure.  For  example,  among  the  Group-VIIB  elements  of 
the  periodic  table,  the  first-row  transition  metal,  iron,1n  its  bcc  crystalline 
form  (a-Fe)  is  exclusively  ferromagnetic,  whereas  the  second-and  third-row 
metals;  ruthenuim  and  osmium,  which  have  the  hep  structure,  are  exclusively 
superconducting.  This  can  be  understood  directly  and  simply  in  terms  of  the 
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diametrical  cluster  molecular-orbital  criteria  for  these  two  phenomeon.  As 
already  discussed  in  Section  III,  for  superconductivity  to  occur,  there  must 
be  degenerate  or  nearly  degenerate  molecular  orbitals  at  the  Fermi  energy 
with  spatially  extended  bonding  components,  e.g.  the  "layered"  pit-bonding 
components  of  A1  and  (TMTSF^PFg  shown  in  Figs.  6  and  25,  respectively,  and  the 
"tubular"  d6-bonding  components  shown  in  Figs.  11  through  14.  In  contrast, 
for  ferromagnetism  or  local  magnetic  moment  formation,  as  already  pointed  out 
in  Section  I,  there  must  be  spatially  localized,  exclusively  anti  bonding  molecular 
orbitals  at  the  Fermi  energy  in  the  non-spin-polarized  limit,  as  exemplified 
for  iron  in  Fig.  14  of  Ref.  41  and  for  Mn  impurities  in  aluminum  in  Figs. 

9  and  10.  Indeed,  it  was  shown  in  Section  III  how  one  can  understand 

the  destructive  effect  of  such  magnetic  impurities  on  the  superconductivity  of 
pure  aluminum  simply  in  terms  of  the  emptying  of  the  aluminum  pn-bonding 
molecular  orbitals  at  the  Fermi  energy  (Figs.  4  through  6)  and  the  occupation 
of  purely  antibonding  Mn-Al  dpc*  and  dp^*  orbitals  [Figs.  9  and  10], 
which  are  responsible  for  magnetic  spin  fluctuations.  Recent  SCF-Xoc-SW 
molecular-orbital  studies  of  the  electronic  structures  and  magnetism  of  iron 
clusters  [41]  indicate  that  a  high  density  of  localized  antibonding  molecular 
orbitals  at  the  Fermi  energy  in  the  non-spin-polarized  limit  is  consistent  with 
the  Stoner  criteria  for  the  occurence  of  ferromagnetism,  namely,  that 
a  paramagnetic  metal  such  as  iron  will  be  unstable  with  respect  to  ferromagnetic 
ordering  if 

I(eF)N(eF)  >  1.  (65) 

where  N(eF)  is  the  electronic  density  of  states  at  the  Fermi  energy  of  the 
paramagnetic  metal  and  I(eF)  is  the  Stoner  exchange  parameter.  In  Fig.  15, 
it.  is  clear  that  the  Fermi  energy  of  paramagnetic  iron  is  coincident  with  a 
region  of  high  density  of  states  corresponding  to  antibonding  cluster  molecular 
orbitals.  In  Ref.  41,  it  is  shown  that  these  antibonding  orbitals  also  lead 
to  large  values  of  I(cF)  because  of  the  spatial  localization  of  the  wave 


-56- 


functions  in  regions  between  the  antibonding  nodes  and  atomic  nuclei 
(see  Fig.  14  of  Ref.  41 ) . 

The  case  for  anti  ferromagnetism  vs.  superconductivity  deserves  separate 
discussion.  It  has  been  shown  in  Section  III  [see  Eqs.  (6)  and  (7)]  that 
the  "valence-bond-like"  singlet  configuration  of  "layered"  pu-bonding.  or 
"tubular"  d6-bonding  molecular-orbital  components  and  4»_  at  the  Fermi  energy 
in  a  superconductor,  which  are  do*  or  dir*  antibonding  between  second-nearest- 
neighbor  atoms  along  directions  perpendicular  to  the  layers  or  tubes,  corresponds 
to  an  "antiferromagnetic"  ordering  of  spatially  delocalized  spins  on  the 
alternating  layers  or  tubes.  This  spin  configuration  is  somewhat  similar  to  the 
conduction-band  collective  anti  ferromagnetism  or  spin-density  wave  predicted  by 
0verhauser[90]  to  occur  in  certain  metals  and  to  be  lower  in  energy  than 
the  paramagnetic  state.  Indeed,  there  is  recent  experimental  evidence  for  an 
"anti ferromagnetic"  ordering  of  conduction  electrons  in  A15  superconductors  [80] 
(see  also  Section  IVC).  On  the  other  hand,  recent  SCF-Xa-SW  cluster  models 
for  the  electronic  structure  of  chromium,  [111]  including  the  effects  of 
magnetic  spin  polarization,  suggest  that  conventional  anti  ferromagnetism  in 
this  material  is  associated  with  molecular  orbitals  at  the  Fermi  energy  which 
have  spatially  localized  d-r*-antibonding  character  between  nearest-neighbor 
atoms  and  spatially  delocalized  dc-bonding  character  between  second-nearest 
neighbors  (see  Fig.  20).  Thus  while  both  Cooper-Schafroth  pair  formation 
and  conventional  anti  ferromagnetism  are  associated  with  composite  bonding¬ 
antibonding  molecular-orbital  topology  at  the  Fermi  energy,  the  two  phenomena 
are  distinguished,  at  least  for  the  systems  investigated  thus  far,  by  different 
spatial  distributions  of  this  composite  character  between  the  atoms. 

Thus  the  SCF-Xa-SW  cluster  molecular-orbital  approach  can  be  considered 
to  provide  a  simple  unifying  chemical  picture  for  understanding  the  incidence 
of  superconductivity,  ferromagnetism,  and  anti  ferromagnetism  entirely  on  the 
basis  of  the  molecular-orbital  topology  at  the  Fermi  energy. 
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VII.  METALS  THAT  ARE  NEITHER  SUPERCONDUCTING  NOR  MAGNETIC 

The  molecular-orbital  criteria  for  superconductivity  and  magnetism  des¬ 
cribed  in  the  preceding  sections  can  also  be  employed  to  explain  some  metals 
are  neither  superconducting  nor  magnetic.  Among  these,  the  noble  metals  Cu, 
Ag,  and  Au  have  been  the  subject  of  much  discussion  concerning  the  possibility 
of  observing  superconductivity  at  extremely  Tow  temperatures.  [112]  SCF-Xa-SW 
studies  of  a  19-atom  copper  cluster  representing  bulk  crystalline  copper 
[39]  suggest  that  superconductivity  is  unlikely  to  occur  at  any  temperature 
(or  reasonable  pressure).  The  Fermi  energy  coincides  with  a  molecular 
orbital  of  a^  "spherical1’  symmetry  which  is  primarily  of  the  so* -anti bonding 
type  and  which  therefore  lacks  the  coherent  extended  spatial  character 
necessary  for  superconductivity  (see  Fig.  4  of  Ref.  39  for  a  picture  of  this 
orbital).  Uhi’e  there  is  some  antibonding  d-orbital  contribution  to  this 
orbital  due  to  the  hybridization  effect  of  the  filled  d  band  below  the 
Fermi  energy,  it  is  too  small  and  diffuse  to  give  rise  to  localized  or 
collective  magnetic  moments  in  pure  copper.  The  only  orbital  that  is 
capable,  in  principle,  of  forming  a  superconducting  state  in  copper  is  an 
unoccupied  t^  p--pc*  orbital  of  the  type  shown  above  (Figs.  4-7)  to  be 
occupied  at  the  Fermi  energy  in  aluminum  and  responsible  for  its  super¬ 
conducting  state.  Thus  the  question  of  superconductivity  in  copper  depends 
on  whether  this  type  of  orbital,  which  lies  well  above  the  Fermi  energy,  can 
be  brought  into  coincidence  with  the  Fermi  energy  by  any  practically  attain¬ 
able  pressure.  Similar  arguments  are  applicable  to  silver  and  gold. 

Another  interesting  example  where  there  has  been  much  speculation  about 
the  possibility  of  superconductivity  is  metallic  hydrogen  at  very  high 
pressure.  [113-115]  SCF-Xa-SW  molecular-orbital  studies  [98]  suggest  that, 
while  the  electronic  structure  of  solid  atomic  hydrogen  at  high-pressure 
Interatomic  distances  and  crystal  structure  is  consistent  with  metallic 
behavior,  this  substance  is  unlikely  to  be  superconductor.  The  Fermi  energy 


coincides  with  an  so*  antibonding  molecular  orbital  of  the  type  shown  above 
to  be  associated  with  the  nonsuperconducting,  nonmagnetic  state  of  copper. 
Thus,  as  in  the  example  of  copper,  the  occurrence  of  superconductivity  in 
metallic  hydrogen  depends  on  whether  an  unoccupied  spatially  extended 
prr-bonding  orbital,  which  lies  well  above  the  Fermi  energy,  can  be  brought 
into  coincidence  with  the  Fermi  energy  at  any  attainable  pressure. 
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Fig.  1 

Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


Fig.  6 


Fig.  7 


Fig.  8 


FIGURE  CAPTIONS 

Schematic  representation  of  the  "real  space"  molecular  orbitals  and 
corresponding  "t-space"  band  structure  of  a  benzene  molecule. 

Comparison  of  the  SCF-Xa-SW  molecular-orbital  energy  levels  of  a 
43-atom  aluminum  cluster,  representing  the  local  environment  of 
fee  crystsalline  aluminum  up  to  ttvird-nearest  neighbors,  with  the 
bulk  energy  bands. 

Comparison  of  the  SCF-Xa-SW  molecular-orbital  energy  levels  of  a 
43-atom  aluminum  cluster,  representing  the  local  environment  of 
crystalline  aluminum  up  to  third-nearest  neighbors,  with  the  measured 
X-ray  photoelectron  spectrum  (XPS)  of  an  aluminum  crystal. 

Contour  map  of  the  wave  function  for  the  highest  occupied  molecular 
orbital,  tj  (sp),  a  43-atom  aluminum  cluster,  plotted  in  the  (200) 
crystallographic  plane  up  to  second-nearest  neighbors.  The  solid 
and  dashed  contours  represent  positive  and  negative  values,  respectively, 
of  the  wave  function. 

Contour  map  of  the  wave  function  for  the  highest  occupied  molecular 
orbital,  t^(sp),  of  a  43-atom  aluminum  cluster,  plotted  in  the  (110) 
crystallographic  plane  up  to  second-nearest  neighbors. 

Three-dimensional  contour  map  of  the  wave  function  for  the  highest 
occupied  molecular  orbital,  t]u(£p)>  °f  a  43-atom  aluminum  cluster,  i.e. 
the  three-dimensional  version  of  the  contour  maps  of  Figs.  4  and  5. 
Schematic  representation  of  the  spatially  extended  wave  function  of 
the  highest  occupied  molecular  orbital,  t-^Up),  of  an  aluminum  cluster, 
plotted  in  the  (200)  crystallographic  plane,  when  the  cluster  size  is 
increased  to  include  many  atoms  in  the  crystal. 

SCF-Xa-SW  molecular-orbital  energy  levels  of  a  43-atom  cluster 
representing  the  local  environment  of  a  Mn  impurity  in  an  aluminum 
crystal  up  to  third-nearest  neighbors.  The  calculated  levels  are 


compared  with  the  measured  X-ray  photoelectron  spectrum  (XPS)  of  bulk 
aluminum  containing  a  dilute  concentration  of  substitutional  Mn 
impurities. 

Fig.  9.  Contour  map  of  the  highest  occupied  cluster  molecular-orbital 

wave-function  for  a  Mn  impurity  in  an  aluminum  crystal. 

Fig. 10.  Contour  map  of  the  highest  occupied  eg  cluster  molecular-orbital  wave 

function  for  a  Mn  impurity  in  an  aluminum  crystal. 

Fig. 11.  Schematic  perspective  drawing  of  "tubular"  dd-bonding  molecular- 

orbital  components  along  a  line  of  atoms  arising  from  the  overlap 

of  d  atomic  orbitals  along  the  z-axis. 
xy 

Fig. 12.  Schematic  perspective  drawing  of  "tubular"  dd-bonding  molecular- 

orbital  components  along  a  line  of  atoms  arising  from  the  overlap 

of  d  2  2  atomic  orbitals  along  the  z-axis. 
x  -y 

Fig.lJ.  Schematic  perspective  drawing  of  "tubular"  dd-bonding  cluster 

molecular-orbital  components  in  a  bcc  crystal  arising  from  the  overlap 

of  d  atomic  orbitals, 
xy 

Fig. 14.  Schematic  perspective  drawing  of  "tubular"  dd-bonding  cluster  molecula 

orbital  components  in  a  bcc  crystal  arising  from  the  overlap  of 

d  2  2  atomic  orbitals, 
x  -y 

Fig. 15.  Composite  molecular-orbi tal  density-of-states  profile  for  bcc 
transition  metals. 

Fig. 16.  Schematic  perspective  drawing  in  a  plane  of  three  atoms  of  "layered" 
bonding  molecular-orbital  components  arising  from  the  overlap  of 
d22  atomic  orbital  annular  regions  (the  dashed  contours). 

Fig. 17.  Contour  map  of  the  wave  function  for  the  highest  occupied  tjg 
(d  ,  d  ,  d  )  cluster  molecular  orbital  of  vanadium 

xy  yj  Xfc 
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Fig. 13. 

Fig. 19. 

Fig. 20. 

Fig. 21. 

Fig. 22. 

Fig. 23. 

Fig. 24. 

Fig. 25. 

Fig. 26 

Fig.27(a) 

Fig. 27(b) 


Contour  maps  in  two  different  planes  of  the  wave  function  for  the 

highest  occupied  e  (d  2  2)  cluster  molecular  orbital  of  niobium. 

9  *  -y 

Contour  map  of  the  wave  function  for  highest  occupied  eg (dz2 )  cluster 
molecular  orbital  of  niobium. 

Contour  map  of  wav§  function  for  the  highest  occupied  eg(dz2)  cluster 
molecular  orbital  of  chromium. 

Contour  map  of  the  wave  function  for  the  highest  occupied  cluster 
molecular  orbital  of  zirconium. 

Contour  map  of  the  wave  function  for  the  lowest  unoccupied  cluster 
molecular  orbital  of  zirconium. 

Contour  map  of  the  wave  function  for  the  d6-bonding  molecular  orbital 
at  the  Fermi  energy  of  Nb^Sn,  shown  along  four  atoms  of  a  Nb  chain. 
Schematic  representation  of  the  wave  function  for  the  molecular 
orbital  at  the  Fermi  energy  of  PdHx  (x  ^  1). 

Schematic  representation  of  the  wave  function  for  the  molecular 
orbital  at  the  Fermi  energy  of  (TMTSF^PFg  in  relation  to  the  crystal 
structure. 

Crystal  structure  of  TTF-TCNO. 

Contour  map  of  the  wave  function  for  the  cluster  molecular  orbital 
at  the  Fermi  energy  of  amorphous  Cu-Zr. 

Contour  map  of  the  wave  function  for  the  cluster  molecular  orbital 


at  the  Fermi  energy  of  amorphous  Pd-Si . 


FIGURE  2 


-73- 


COUNTS  (ARB.  UNITS) 


FIGURE  3 


BINDING  ENERGY  (ev) 


74 


FIGURE  4 


vy  v 


FIGURE  6 


ENERGY (RYDBERGS) 


FIGURE  10 


i  . 


AD-A119  784  MASSACHUSETTS  INST  OF  TECH  CAMBRIDGE  CENTER  FOR  MATE— ETC  F/G  11/6  \ 

molecular-orbital  basis  for  SUPERCONDUCTIVITY  in  HIGH-  AND  LOW— ETC(U) 

SEP  82  K  H  JOHNSON,  R  P  MESSMER  N00014-81-K-0499  ' 

UNCLASSIFIED  TR-4  NL 


I 


FIGURE  23 


